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ABSTRACT 


\ 


his  dissertation  presents  the  quantitative  determination  of  the 
applicability  of  a likelihood  ratio  processor  in  the  classification 
of  sonar  returns.  The  ability  of  a likelihood  ratio  processor  to 
classify  simple,  similar  targets  was  measured  under  highly  con- 
trolled conditions , The  five  targets  used  were  7 in.  diam,  5 in. 
diam,  and  5 in.  diaA,  solid  aluminum  spheres,  a 5 in.  diam  styrofoam 
sphere,  and  a 5 in.  cU.am  hollow  aluminum  sphere.  The  two  trans- 
missions used  were: 

(l)  a sinusoidal  pulse  with  a frequency  of  70  kHz  and  a 
pulse  length  of  80  psec,  and 
(?)  a linear  PM  of  bandwidth  12.5  kHz  with  a center 

frequency  of  6j.75  kHz  and  a pulse  length  of  3 msec. 


Since  the  S/N  ratios  of  the  recorded  echoes  were  large,  Gaussian 
noise  from  a random  noise  generator  was  added  to  the  echoes  in  order 
to  simulate  a more  realistic  problem  in  target  classification.  In 
classifying  messages  (echo-noise  combinations)  from  two  target  classes, 
the  performance  of  the  likelihood  ratio  processor  was  measured  by 
means  of  ROC  curves.  The  classification  was  based  on  the  quadrature 
components  of  the  messages.  For  the  sinusoidal  pulse,  the  likelihood 
ratio  processor  performed  satisfactorily  in  classifying  the  messages 
for  a majority  of  the  target  pairs  when  the  S/iJ  ratio  of  the  messages 
was  0 dB  or  larger.  ^ ' SC-"" 

For  the  sinusoidal  pulse,  the  effect  of  variations  in  the  S/^  ratio 
of  the  messages  on  the  performance  of  the  likelihood  ratio  processor 
was  measured.  The  S/N  ratio  was  found  to  be  a sensitive  parameter 
that  influenced  the  perfonnance.  Also  for  the  sinusoidal  pulse,  the 
likelihood  ratio  processor  was  able  to  classify  target  classes  using 
only  the  in-phase  or  out-of- phase  components  of  the  messages. 


The  effect  on  performance  of  an  uncertainty  in  epoch  over  a cycle  of 
the  carrier  was  measured  for  the  sinusoidal  pulse.  The  uncertainty 
in  epoch  brought  about  the  same  change  in  performance  as  a 6 dB  drop 
in  the  S/TI  ratio  of  the  messages . 


Since  the  reciprocal  of  the  pulse  length  of  the  sinusoidal  pulse 
was  equal  to  the  bandwidth  of  the  linear  TM,  it  was  asserted  and 
found  that  equivalent  performances  could  be  obtained  with  either 
of  the  two  transmissions  if  the  following  equation  were  true: 
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n '\ 
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where 


E is  the  signal  energy  in  an  echo, 

< n2  > is  the  mean  square  of  the  added  Gaussian  noise, 
sp  denotes  a sinusoidal  pulse,  and 
M denotes  a linear  M. 

E 

The  dependence  of  performance  on  p — implies  that  a linear  M 

< N > 

has  more  potential  for  improvement  in  the  performance  of  a likelihood 
ratio  processor  than  a sinusoidal  pulse,  since  the  energy  in  a 
sinusoidal  pulse  is  restricted  by  the  effects  of  cavitation. 
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CHAPTER  I 


INTRODUCTION 

The  objective  of  this  investigation  is  to  determine 

quantitatively  the  applicability  of  a likelihood  ratio  processor  in 
» 

the  classification  of  sonar  returns  from  targets. 

The  likelihood  ratio  processor  implements  the  Bayes' 
decision  rule  in  the  classification  of  targets.  The  Bayes'  decision 
rule  is  optimum  in  the  sense  that  no  other  rule  will  yield  a smaller 
average  risk  to  the  observer. 

The  ability  of  a likelihood  ratio  processor  to  classify 
correctly  sonar  returns  from  similar,  simple  targets  was  measured 
for  both  a sinusoidal  pulse  and  a linear  M transmission.  The  classi- 
fication of  the  targets  was  based  on  the  in-phase  and  out-of-phase 
quadrature  components  of  sonar  returns . Measurements  were  made  on  the 
ability  of  the  likelihood  ratio  processor  to  classify  correctly  sonar 
returns  using  the  in-phase  components  and  the  out-of-phase  components 
separately.  In  addition,  the  effects  of  variations  in  the  average 
S/N  ratio  and  the  epoch  (or  reference  time)  of  the  sonar  returns  on 
the  ability  of  the  likelihood  ratio  processor  to  classify  correctly 
were  measured. 

A discussion  of  the  statistical  decision  theory  and  the 
application  of  theory  in  psychophysical  experiments  is  presented  in 
Signal  Detection  Theory  and  Psychophysics  by  David  Green  and 
John  Swets . The  techniques  of  specification  of  an  optimal  detection 
process  and  of  the  use  of  ROC  curves  as  a measure  of  experimental 
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observations  are  similar  to  those  implemented  in  target  classification. 
In  the  classification  of  similar,  simple  targets,  a metric  is  applied 
as  a geometrical  description  of  target  echoes . A discussion  of  a 
measurement  of  the  similarity  of  points  in  a multidimensional  vector 
space  expressed  by  a metric  is  found  in  Decision-Making  Processes  in 
Pattern  Recognition  by  George  Sebestyen, 

The  problem  of  target  classification  is  important  in  many 
areas  of  information  processing.  The  use  of  a likelihood  ratio 
processor  in  target  classification  is  applicable  to  sonar  and  radar 
surveillance,  to  geophysical  prospecting,  to  classification  of  magnetic 
spin  echoes,  and  to  radio  astronomy. 
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CHAPTER  II 


THEORY 


A.  Likelihood  Ratio 

The  following  discussion  concerning  likelihood  ratio  assumes 
that  an  active  source  transmits  a bandlimited  signal  directed  at  one 
of  two  possible  targets  present  in  a propagating  medium. 

The  signals  from  the  targets  are  defined  as 

S'^(t;cp^),  j = A or^  , 

where 

j = A refers  to  a signal  from  target  A,  and 
j = B refers  to  a signal  from  target  B. 

The  starting  phase  of  a signal  is  cp^and  it  is  contained  within  the 
interval  -n  < 9^  < « . The  assumed  density  function  for  9^  is 
P (cp^)  = This  assumption  is  made  since  the  epoch  of  a signal  at 

the  receiver  is  not  known.  The  receiver  output  is  defined  as  V(t), 
the  echo, 
where 

V(t)  = N(t)  + (ticp^)  . (2.1) 

The  echo  V(t)  is  bandlimited  to  a bandwidth  W with  a center 
frequency  f^.  The  additive  noise,  N(t),  is  composed  of  ambient  noise 
in  the  medium  and  in  the  receiver;  l.e.,  reverberation  is  not  considered. 


Since  the  echo  V(t)  is  an  analog  signal,  it  must  be  digitized 
by  an  analog-to-digital  converter  before  it  is  used  in  any  digital 
computations.  After  V(t)  is  sampled  and  digitized,  it  is  expressed 
as  = (v(iAt))  = (V^).  The  time  interval  between  the  samples  is  At, 
and  I = 1,  2,  ....  TF,  where  TF  is  the  total  number  of  samples  of  V(t) 
that  are  sufficient  to  reconstruct  V(t)  from  a specified  interpolation 
formula.  F is  the  sampling  rate  used  to  obtain  samples  of  V(t),  and 
T is  the  pulse  length  of  V(t).  The  sampled  vector  ^ defines  a point 
in  a TF- dimensional  observation  space  F. 


In  this  study  the  noise  N is  restricted  to  a stationazy 
irocess  of  zero  mean.  Also,  I?  is  bandlimited  to  a 


bandwidth  W with  center  frequency  f . The  multivariate  probability 


density  function  describing  TF  discrete  samples  of  the  noise  is 
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probability  of  and  an  a priori  probability  of  where 

= 1.  If  the  a priori  probabilities  q^  and  q^  are  not  known, 
It  Is  necessary  to  train  the  observer  with  a representative  sample 
set  in  order  to  obtain  values  for  q^  and  q^.  The  hypothesis  is 
tested  by  observing  an  echo  After  observing  is  either 

accepted  or  rejected,  A decision  rule,  which  is  needed  to  make  one 
of  the  two  possible  decisions,  will  be  formulated  by  dividing  the 
observation  space  F into  two  disjoint  subspaces  and  F^,  where 
F = F^UFg.  A decision  is  made  to  accept  if  ^ is  contained  in  F^, 
or  reject  if  ^ falls  into  Fg. 


F = F^UFg 


F^nFg  = 0 


Since  many  choices  for  the  boundary  are  possible,  a method 
is  needed  for  determining  an  optimum  decision  boundary.  In  determining 
an  optimum  decision  boundary,  it  is  possible  to  make  two  kinds  of  errors; 

(1)  rejection  of  when  it  is  true,  and 

(2)  acceptance  of  when  it  is  false. 
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A choice  of  a bouridaiy  determines  the  probability  of  making  an  error 
of  each  type  as  follows: 


K*?  ) " /*  ~ probability  of  error  No.  1 

° Jr  ° function  of  cp  , where  cp  i 


1 as  a 

is  the  (2.6) 


starting  phase,  and 


-1/2  a^  (y  -S^(cp  )')  (v  -S^(cp  ) 
' mn\  m m'^oV  \ n n'  o 


■ -^2 


m,n=l 


) - / P("^|Hp;cp  )d^  = probability  of  error  No.  2 as  a 

^ r ^ fnnr»+.-t  rin  rtf'  m 


function  of  cp  , 
o 


where 


- -1/2  V*  ^ -S®(cp  (y  -S®(cp  ) 

1/2  ' 2^  mn\  m m'^oV  \ n n'^o  , 


m,  n=l 


It  is  possible  to  make  either  Ct{(p^)  or  6(9^)  arbitrarily 

small  by  adjusting  the  boundary  to  decrease  the  number  of  elements 

in  either  or  F^,  but  it  is  not  possible  to  make  Q!(<Pq)  s-nd  6(9^) 

2 

both  arbitrarily  small. 

A general  solution  to  the  hypothesis  testing  problem  known 
as  the  Bayes'  solution  is  possible  if  a set  of  cost  functions 
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which  describe  the  consequences  of  deciding  that  hypothesis  is 
true  when  H.  is  true  (i, j=l,2),  is  available.  In  this  study,  the 
four  cost  functions  are  predetermined  as  follows: 

= 0 = cost  to  the  observer  if  he  decides 
is  true  when  it  is  true, 

^12  ~ cost  to  the  observer  if  he  decides 
is  true  when  is  true, 

0^2^  = cost  to  the  observer  if  he  decides 
is  true  when  is  true, 

C22  = 0 = cost  to  the  observer  if  he  decides 
H2  is  true  when  it  is  true. 

The  Bayes'  decision  rule,  which  results  in  the  minimum 
average  risk  to  the  observer  and  establishes  a boundary  between 
and  Fg,  is  formulated  by  defining  two  conditional  expected  risk 


functions  as  follows: 


~ expected  risk  given  is  true 


~ expected  risk  given  Hg  is  true 


The  average  risk  is  given  by 


I / 


f 

I 


4. 

t * 


R = 


/j: 


P(^lH^jcp^)  P(cpQ)dcp^d^ 


•''“/Jj 


From  the  result  that 


-U  P(^|H^;q)^)  P(cp^)dcp^d^  , or 
1 = j P(cp^)cLcp^d^ 

• II 


P(^lH^jcp^)  P(cp^)dcp^d^  ^ 


Eq.  (2.13)  becomes 


(2.L3) 


P(^|H2;<P^)  P(cp^)d(p^d^  . 


(2.14) 


R 'li^21 


•'^l  (2.15) 

Since  the  first  term  in  Eq.  (2.I5)  is  a constant,  the 
average  risk  is  a minimum  when  the  integral  over  is  a minimum. 

Since  the  integrand  consists  of  two  non-negative  terms,  the  integral 
attains  its  minimum  value  if  the  following  two  conditions  are  satisfied: 
1)  The  subspace  F^  consists  of  those  values  of  ^ for  which 


’Ai  rp(v’IHi;<p„)  > <l2Cl2  f 

J -Tt  J -n 


(2.16) 


^ 1 hi 


J 
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2)  The  subspace  consists  of  those  values  of  ^ for  which 

rn  rn 

V2I  / P('Po)'i'Po  < <12^12  / ^(^iHpJcp^)  P(<PQ)d'Po  • (2.17) 


The  generalized  likelihood  ratio^,  A(^,  is  defined  as 


/It 

P(^|H2;cp^)  P('Po)dcp^ 

2 

rst 

J -7t 


(2.18) 


In  summary,  the  Bayes'  decision  rule  assigns  to  P^  those 

j j ^1^21 

V for  which  A(Vj  < — - — , and  assigns  to  P those  V for  which 

V12  ^ 

J ^1^21  *^1^21 

^ — r — • number — ^ — is  interpreted  as  a threshold.  A 

‘^2^12  V12 

device  that  implements  the  generalized  likelihood  ratio  calculations 
in  hypothesis  testing  is  referred  to  as  a generalized  likelihood  ratio 


processor. 

If  the  epoch  of  the  echoes  were  known,  the  probability 
density  function  for  the  starting  phase  would  become  P(9^)  = &(<Pq)> 
a Dirac  delta  function,  and  the  generalized  likelihood  ratio  would 
reduce  to  the  likelihood  ratio  , which  is  defined  as 


i(^  = 


P(^|H2) 

p(^|h^) 


^Ibld.,  p.  58. 


Ibid.,  p.  59. 
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Similarly,  a device  that  implements  likelihood  ratio 
calculations  in  hypothesis  testing  is  called  a likelihood  ratio 
processor.  The  threshold  used  in  likelihood  ratio  calculations  is 

%^21 

^2^12  ■ 

Upon  the  substitution  of  Eq.  (2.7)  and  Eq.  (2.9)  into 
Eq.  (2,l8),  the  generalized  likelihood  ratio  becomes 


, -1/2  y (v  - S^(cp))  (v  - S®(cp) 

1/2  ^ mn  \ m V \ n n'^  , 


-1  a 

^ I ' mn' 

2it 


/ 


m,n=l 


dcp 


TF 


hf 


1/2  E 


p 

( V 

mn  \ m m 


mn ' 


m,n=l 


{2n) 


TF/2 


dcp 


Upon  the  simplification  of  Eq.  (2.20),  a(^  becomes 


a(7) 


f 


y fv 


(cp)  - 1/2  s^(cp)s; 


m,n=l 


?(.)] 


dcp 


TF  r -I 

a^  V S^(cp)  - 1/2  S^(cp)S^(cp) 
7 . mnl  m n'^^  ' vy  ' n'^^ 

m,n=l  ^ 


dg 


(2.20) 


(2.21) 


I; 

i I 


f 

p 

r '» 


, t.  Mi,*. 


Upon  the  substitution  of  Eq.  (2.7)  and  Eq.  (2.9)  into 


Eq.  (2.19)  and  simplification,  the  likelihood  ratio  i(V)  for  cp^  = 0 


becomes 


Y'  [v  S®  - 1/2  S®S®1 

/ . mn  m n ' m n I 


1/2 


m,n=l 


TP  r T 

a^  V - 1/2 
/ . mn  I m n ' m n 


(2.22) 


m,  n=l 


where 


S®(t  ;cp  = 0)  = S^(t  ) = S 


, and 


S^(t  ;cp  = 0)  = S^(t  ) = 

'■m  o 'nr'  m 


m' 


m 


B.  Narrow  Band  Signals 


Many  signals  used  in  physical  applications  are  bandlimited 


to  a frequency  interval  I f 


W 


< f < f + 


where  W is  small 


compared  to  f^.  Signals  of  this  type  are  referred  to  as  narrow  band 


signals.  A representation  for  these  narrow  band  signals  is 

S(t)  = ST(t)  cos  2jtf  t + S_(t)  sin  2jtf  t 
' ' 1'  o 2.^  ’ o 


(2.23 


where 


S^(t)  = A(t)  cos  cp(t)  , and 
S2('t)  = -■A(t)  sin  cp(t) 

Sji^(t)  and  S2(t)  are  bandlimited  to  the  frequency  interval  | f | < 

A(t)  is  the  envelope, 

cp(t)  is  the  phase,  and 

f^  is  the  carrier  frequency  of  S(t). 

For  the  present  ST(t)  and  S..,(t)  are  assumed  to  be  of  infinite  duration 


r 


S^(t)  and  82(0)  are  called  the  quadrature  components  of  the 
signal  S(t),  These  functions  contain  both  the  envelope  and  phase 
information  of  S(t),  and  are  of  lower  frequency  than  the  signal  S(t). 
S^(t)  is  referred  to  as  the  in-phase  component  and  S2(t)  is  referred 
to  as  the  out- of- phase  component. 

It  is  assumed  that  S(t),  S^(t),  and  SgCt)  have  Fourier 
transforms.  The  Fourier  transforms  of  S(t),  S^(t),  and  S2(t)  are 
denoted  as  S(f),  S^(f),  and  32(f).  Graphically,  S(f),  S^(f),  and  32(f) 
can  be  depicted  in  the  frequency  domain  as  shown  in  Fig.  2.1. 

The  frequency  domain  occupied  by  S^(f)  and  32(f)  is  con- 
structed by  heterodyning  3(f)  to  a center  frequency  of  0. 

The  function  3^(f)  can  be  represented  over  the  frequency 
range  - W/2  < f < w/2  by  the  complex  Fourier  series 


i2irnf 


yf)  . 2:  0^  e “ 


(2.24) 


where 


i2Ttnf 


: = ■=  / s,  (f)  e ^ 

W^-W/2  1 


(2.25) 


n = 0,  ±1,  ±2,  ..., 


».  Mit  k.  . •t: 


f 

I 


t; 

i' 

• 


r i 
t 
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/\ 

The  infinite  set  of  coefficients  determine  S^(f)  uniquely. 
The  Fourier  transform  of  S^(f)  can  be  written  as 


Si(t)  = 

^ '-W/2 


/■w/2 


i2nft 


df 


(2.26) 


Comparison  of  Eq.  (2,25)  and  Eq.  (2.26)  yields 

■ I ®i(S) 


(2.27) 


Therefore  the  in-phase  component  S^(t)  and  the  out- of- phase 
component  S2(t)  must  be  sampled  at  a minimum  rate  of  W 1'^  order 


to  extract  all  of  the  information  contained  in  S^(t)  and  S2(t). 


Samples  of  S,  (t)  obtained  at  the  minimum  rate  of  W will  be 

1'  sec 


independent  of  each  other.  Namely,  a change  in  one  or  more  of  the 


sampled  values  S^(n/w),  n = 0,  ±1,  ±2,  ...,  does  not  alter  the  rest,' 


Samples  of  Sg(t)  obtained  at  the  minimum  rate  of  W also  will 


be  independent  of  each  other. 

It  was  assumed  that  S(t)  was  of  infinite  duration.  However, 
if  the  duration  of  a signal  is  not  infinite,  but  is  long  compared  to 
its  fluctuation  time  (l/2w),  an  interpolation  fomula  can  represent 


the  signal  from  its  samples  in  the  interval  (-T/2,  T/2)  to  within 

6 


an  error  that  is  of  the  order  l/l . 

If  the  pulse  length  T of  the  quadrature  components  is 
finite,  a minimum  of  TW  samples  of  each  quadrature  component  is  re- 
quired for  digital  processing  of  the  quadrature  components. 


■^David  Middleton,  An  Introduction  To  Statistical  Communication  Theory, 
McGraw-Hill  Book  Co.,  New  York,  N.  Y.,  I96O,  p.  209. 


Ibid.,  p.  211. 


4 
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C . Quadrature  Sampling 

7 

The  technique  of  quadrature  sampling  is  used  to  obtain  the 
quadrature  components  of  a narrow  band  signal.  The  use  of  quadrature 
sampling  to  obtain  the  quadrature  components  eliminates  the  need  to 
heterodyne  a signal. 

The  in-phase  component  S^(t)  of  S(t)  is  obtained  by  sampling 

S(t)  at  times  t = nK/f^,  n = 1^  2,  ...  TW^.  TW^^  is  the  number  of  samples 

of  S (t)  where  W = f /K.  The  sampling  parameter,  K,  is  a positive 
1.  o 

integer  that  determines  the  time  interval  between  samples  of  S^(t)  in 
multiples  of  sampling  parameter,  K,  is  equal  to  [f^/w] 

where  [f^/w]  denotes  the  largest  integer  K S 

To  obtain  the  out-of -phase  components,  S^Ct),  S(t)  is  sampled 
at  times  t = ~ '^>  ‘ ‘ ‘ ^ small  compared 

to  f^,  the  approximation  S^CnK/'f^  + s S^CnK/f^)  can  be  made. 

Otherwise  S^CnK/f^)  may  be  computed  from  S2(nK/f^  + l/lf^)  by  using 
an  interpolation  formul.a. 


From  the  samples  of  S^(t)  and  S2(t),  the  envelope  and  phase 
of  S(t)  can  be  computed  at  a time  >^K/f^  as  follows: 


/nK  \ _ 


L2fnK  ^ ^ .2fnK 


, and 


(2.28) 


iiiV  t 

2VV 

/nK\ 


(2.29) 


0.  D.  Grace  and  S,  P.  Pitt,  "Quadrature  Sarnpling  of  High  Frequency 
Waveforms,"  J.  Acoust.  Soc.  Am.  (1968)  pp,  l453-i454. 
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Each  sample  of  the  in-phase  component  will  be  independent 

of  other  samples  of  the  in-phase  component  when  = W.  Each  sample 

H 

of  the  out-of-phase  component  will  be  independent  of  other  samples  of 
the  out-of-phase  component  for  this  same  condition.  Also  the  quadra- 
ture components  resolved  from  a signal  are  dependent  on  the  choice  of 
origin  for  the  time  scale. 

If  samples  of  S(t)  cannot  be  obtained  at  the  times 
nK/ f^,  physical  limitations 

on  the  analog-to-digital  converter,  but  can  be  obtained  at  times 


nK/f^,  nK/f^  + n = 1,  2,  ...  TW^,  it  is  permissible  to  redefine 

the  center  frequency  (f^)  as  f^  and  the  bandwidth  (W)  as  W . The 
redefined  center  frequency  f'  is  as  near  the  value  of  f^  as 
the  limitations  on  the  analog-to-digital  converter  allow, 

W = W + 2 |f^-f^l  and  = f^/[f^/W'].  The  equation  for  S(t)  then 
becomes 

S(t)  = S'(t)  cos  cn’t  + S;:.(t)  sin  co’t 
1 o a o 

where  o)'  = 2jtf  . If  W is  small  compared  to  f and  If  - f'l  « f . 

0 0 o ' o o o' 

the  approximation  Sg(nK/f^  + = S^(nK/f^)  can  be  made. 

The  maximum  error  (e^)  made  in  the  approximation 

O 

S^CnK/f^  + 1/Uf  ) s Sp(nK/f  ) is  proportional  to  W/f  . The 


o 

Athanasious  PapouLis,  Limits  on  Band-limited  Signals,"  Proc.  IEEE 
No.  10,  (October  196?),  p.  l680. 


l8 


increase  in  error  (Ae  ),  resulting  from  redefining  the  bandwidth 

as  W and  the  center  frequency  as  f^,  is  proportional  to 

(w'/f  - W 'f  ).  And 
' o o 


v;'  w 


o 


= "a 


f w' 
o 

f 'W 
. o 


The  echoes  discussed  in  Chapter  IV  were  bandlimited  to  a 
bandwidth  of  20  kHz  with  a center  frequency  of  70  kHz.  The  echoes 
were  digitized  at  a rate  of  JOO  kHz  instead  of  the  desired  rate  of 
280  kHz,  Physical  limitations  on  the  analog-to-digital  converter 
necessitated  the  use  of  this  rate.  The  bandwidth  was  redefined  as 
30  kHz  with  a center  frequency  of  75  kHz . The  increase  in  error 
(Ae.)  resulting  from  redefining  the  bandwidth  and  the  center  fre- 
quency  was  0.4  e^. 

In  Fig.  2,2,  the  in-phase  component  S|(t)  and  the  out-of-phase 
component  S2(t)  of  a narrow  band  echo  are  plotted. 

In  Fig.  2.3,  the  envelope  computed  from  the  quadrature 
components  Sj^(t)  and  S2(t)  of  a narrow  band  echo  is  plotted.  The 


ECHO  FROM  5 in.  diom  SOLID  ALUMINUM  SPHERE 


OUT-OF-PHASE  COMPONENT 


FIGURE  2.2 

ECHO  FROM  5 in.  diam  SOLID  ALUMINUM  SPHERE 
AND  ITS  IN-PHASE  AND  OUT-OF-PHASE  COMPONENTS 
= 75  kHz  W'  = 30  kHz 

TRANSMITTING  WAVEFORM  WAS  A SINUSOIDAL  PULSE 
f = 70  kHz  T = 80  /xsec  PEAK  |AMPLITUDE|  = 1.0  in. 


20 
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total  error  (l.^  e ) made  in  the  approximation 

Pl 


is  considered  to  he  small  since  the  envelope  plotted  in  Fig.  2.J 
would  closely  correspond  to  a plot  of  the  actual  envelope  of  the  echo 
if  the  actual  envelope  were  plotted. 

D.  Implementation  of  Livelihood  Ratio  Calculations  Using  (otiuadrature 
Components  Obtained  by  Quadrature  Sampling 

The  quadrature  components  obtained  by  quadrature  sampling  can 

be  used  in  a generalized  likelihood  ratio  calculation.  Pertinent  to 

the  justification  of  the  use  of  quadrature  sampling  is  a generalized 

cross correlation,  such  as 


a(t-u)  V(u)  S'°(t;cp^)  dudt 


The  functions  S'°(t;cp^)  and  V(u)  are  narrow  band  signals  of  bandwidth 
W with  center  frequency  f^,  where  W is  small  compared  to  f^,  and 
a(t-u)  is  a kernel  that  is  a symmetric  function  of  t-u.  The  result 
of  integrating  over  u is  a crosscorrelation 


X(t)  S^(t;cp^)  dt 


where 


X(t)  = I a(t-u)  V(u)  du 


(2.50) 
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Since  X(t)  is  equal  to  a convolution  of  two  functions  in  the  time  domain, 
X(cu)  can  be  described  as 

X(oi)  = a(o3)  V(cd)  (2.51) 

in  the  frequency  domain, 
where 

/\ 

X(cu)  is  the  Fourier  transfom  of  X(t),  etc. 

Since  V(cu)  is  narrow  band,  it  follows  that  X((jd)  is  narrow  band. 

The  aim  of  the  following  arguments  is  to  show  that  the 
generalized  crosscorrelation  in  integral  form  is  approximately  equal  to 
each  of  the  two  following  operations: 


(1) 


where 


F 

a 

mn 


V 

m 


is  a sample  of  the  kernel  a(t-u)  at  the  times  t^  = m/F, 

u = n/F,  and 
n ' ’ 

is  a sample  of  V(t)  at  the  time  t^  = m/F,  etc.,  and 


(2) 


2W^  ^ 

Q P,q=l 


w 

a ^ [V,  S®  (cp) 
pq  Ip  Iq'^"^ 


where 


and 


■^j  is  the  largest  Integer  K = 


Mtti  *i 
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Wn 

a is  a sample  of  the  kernel  a(t-u)  at  the  times 
pq  ^ ^ 


P 


> 

q 


is  a sample  of  V^(t)  at  the  time 


t = 
P 


, etc. 


After  it  has  been  demonstrated  that  the  two  operations  defined 
in  the  last  statement  are  equivalent  to  the  same  generalized  cross- 
correlation,  the  two  operations  can  be  set  equal  to  each  other. 

Thus 


This  demonstration  is  important  since  it  shows  how  the  summations  in 
the  generalized  likelihood  ratio  defined  by  Eq.  (2.2l)  can  be  replaced 
by  summations  involving  quadrature  components . 

A bandpass  signal  such  as  X(t)  is  determined  by  taking  W 
sagp^^s  signal  X(t)  and  W of  the  corresponding  signal 

SGC  SGC 

X(t),  where  X(t)  is  the  Hilbert  transform  of  X(t). 

In  terms  of  samples  of  X(t)  and  of  X(t),  X(t)  may  be 
Q 

written  as 

C.  Van  Schooneveld,  Scsne  Remarks  on  Sampling  Methods  for  a Band 
Pass  Signal,"  Signal  Processing,  NATO  Advanced  Study  Institute, 
(September,  1964),  p.  lj-12. 


2k 


where 


The  functions  b(t--S)  and  d(t-S), 
w w 


are  mutually  orthogonal 


10 


that  is 


, and 


b(t-|)  b(t--3)dt  = ^ Bpq 

<i{t-|)dt  = -^  5pq  , and 
= 0 for  all  p,q 


(2.32) 


(2.35) 


(2.34) 


(2.35) 


The  application  of  the  orthogonal  relations  defined  by 
Eq.  (2.35),  Eq.  (2.34),  and  Eq.  (2.35)  to  the  right  side  of  the 
following  equation, 


^‘^Ibid.,  p.  412. 
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r . /■"  [x(3)  112^  eos  a»f„(t-2) 

1-00  J ~oo 

^ q=-oo 

P=-00 

■*•  ^^(^J'Pq)  Sin  fl(Wt-p)  sin  &tfjt-^)]dt 
3t  ( Wt-  p ) 

results  in  Eq.  (2.37)> 

X(t)S®(t;,pJdt  . ji  ^ [x,S»(.p„)  ^ X,5®(<pJ_ 


(2.56) 


/: 


q=-o 


where 


(2.37) 


= ^q  ’ 

= ^p('Po)  ’ 

Thus  the  crosscorrelation,  / X(t)S^(t;cp  )dt,  is 


ion,  J_^ 


equivalent  to  an  operation  on  the  evenly  spaced  samples  of  X(t), 
X(t)  and  S®(t;cp^),  S®(t;cp^)  obtained  at  the  minimum  rate  W " • 

In  Chapter  II,  Sec.  B,  X(t)  is  written  as 
X(t)  = A(t)  cos  (cp(t)  + oi^t)  = X^(t)  cos  cD^t  + Xg(t)  sin  cc^t 
where 

Xj^(t)  = A(t)  cos  cp(t)  , and 
Xg(t)  = -A(t)  sin  cp(t) 

X^(t)  is  called  the  in-phase  component  of  X(t)  and  X^(t)  is  called 
the  out-of- phase  component  of  X(t). 


t 4. 


'iS*' 
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The  Hilbert  transform  of  X(t)  can  be  obtained  by  shifting 
the  phase  cp(t)  of  X(t)  by  a quarter  of  a cycle  of  the  carrier  f^. 
Thus  X(t)  becomes 


X(t)  = A(t)  cos  ^cp(t)  + u)^t  + -ji  I 

r / / 

= A(t)  cos  q)(t)  cos  I os^t  + j-  sin  cp(t)  sin  + 77!“) 

o'  ^ o'_ 


= X„(t)  cos  o)  t - XT(t)  sin  x-  t 
2'  ^ o 1'  ' o 

If  X(t)  is  sampled  at  the  times  t = nK/f  , n = 0,  ±1, 

no 

±2,  ...  where  K = [f^/w],  the  result  is  samples  of  X^(nK/f^).  If 

X(t)  is  sampled  at  the  times  t = nK/f  + l/^f  , n = 0,  ±1,  t2,  ... 

n o o 

the  result  is  samples  of  X„(nK/f  + 1/^f  ).  If  W is  small  compared 

^'00 

to  then  the  approximation  XgCnK/f^  + 1/^f^)  = X^CnK  f^)  can  be 
made,  and  the  samples  X2(nK/f^),  n = 0,  ±1,  ±2,  ...  can  be  interpreted 
as  samples  of  X(nK'f^),  n = 0,  ±1,  ±2,  .... 

If  samples  of  X^(t),  Xg(t),  and  S^®(t;cp^),  S^ltjcp^)  are 
obtained  at  the  rate  W^,  where  = f ^K,  then  Eq.  (2.57)  can  be 
written  as 
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h 
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In  terms  of  quadrature  components,  Eq.  (2.50)  can  be  written 


as 


X,  (t)  cos  0)  t + X„(t)  sin  ao  t 

I ' ' r\  y ^ ' o 


J a(t-u)  V^(u)  cos  cu^u  + V^(u)  sin  oj^uj  du 


(2.39) 


where 


X(t)  = X^(t)  cos  oD^t  + Xgft)  sin  (a^t  , 


a(t-u)v(u)  = a(t-u) 


V,  (u)  cos  o)  u + V„(u)  sin  cv  u 
1'  o 2'  ' o 


and  X^(t),  X2(t),  a(t-u)V^(u),  and  a(t-u)V2(u)  are  of  bandwidth  W 
with  center  frequency  f^  = 0. 


If  X(t)  is  sampled  at  the  times  t^, 


where 


tg  = qK'f^,  q = 0,  ±1,  ±2,  ...  and  a(t^-u)v(u)  is  sampled 
at  the  times 

Up  = PK/f^,  P = 0,  ±1,  ±2,  . , 


then 


f 

I 


0°  W„ 


Xn(t  ) = ^ V a ^(t  -u  )V,(u  ) , or 

^ p=-00 


<»  w, 


Q 


iq  = 


1 V a '' 

— pq  IP 


W, 


Q p=-“ 


where 


W W W 

Q/4.  ^ Q Q 

a (t-u)  = a =a 

^ q P^  qp  pq 


(2.40) 


(2.41) 


91 


• i /tii  % . . ^ 


A 


If  X(t)  is  sampled  at  the  times  t^  + 

where 

t^  = qK/f^,  q = 0,  +1,  ±2,  ...  and  a(qK/f^  + lAf^,  u)  V(u)  is 
sampled  at  the  times  where  = pK/f^,  p = 0,  ±1,  ±2,  .... 


00  w 

* 5T>  ■ t * fF'  ■ "p  - ir)''2("p  ■" 

^ o P o P o 


If  Wq  is  small  compared  to  f^,  then 

X„(t  + lAf  ) = Xo(t  ) = X_ 

2'  q 'o'  2'  q^  2q  , 

V„(u  + lAf  ) = V„(u  ) = V_ 

2'  p 'o'  2'  p^  2p  , 


and  Eq.  (2.42)  becomes 


1 ” W 

■^2  - -t  E “pq  ''PP  ’ 

p=-oo 


(2.43) 


where 


a ^(qK/f  + l/4f  - pK/f  - l/4f  )=a^(t  -u)=a^=a^ 
''o  'o  ^'o  'o'  ^q  p^  qp  pq 


Substitution  of  Eq.  (2.4l)  and  Eq,  (2.43)  into  Eq.  (2,38) 


results  in 


® W I- 
1 „ Qf-IT 


Q P,q=-« 


X(t)S  (t;cp  )dt  = ~ a ^ S°  (cp  ) + V_  S;;  (cp  ) 

° 2W^  ^ 2p  2q^^o^J  (2,1+k) 
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In  the  previous  equations,  the  sums  on  p and  q are  from 

_oo  to  +“>.  This  assumes  that  the  signals  Vj^(t),  Vg(t),  and  E^(t;rp^), 

So(t;cp  ) are  of  infinite  duration.  Hov/ever,  if  the  duration  of  a 
o 

signal  such  as  V^(t)  is  not  infinite,  hut  is  long  compared  to  its 
fluctuation  time  (l/2V/),  an  interpolation  formula  can  represent  the 
signal  from  its  samples  in  the  interval  (-T/2,  T/2)  to  within  an 
error  that  is  of  the  order 

A representation  of  a signal  such  as  V^(t)  in  the  interval 
( -T/2,  T/2)  by  an  interpolation  formula  is 


V^(t) 


2. 


sin  3t(W_t-q) 
n(W|^t-q) 


(2.45) 


Upon  a translation  of  the  time  coordinate,  namely, 
t = t + T/2  + l/Wj^j  the  interpolation  formula  becomes 


Q / ^ sin  3t(W  t-q) 

' E A t 

q=l  ^ ^ 


(2.46) 


where 


TW,  + 1 = TW 


David  Middleton,  An  Introduction  to  Statistical  Communication  Theor;:) 
McGraw-Hill  Book  Co.,  New  York,  N.  Y. , (19&).  p.  211. 
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Therefore,  if  the  signals  V^(t),  ^ ^t),  and  S^(t;cp^), 
S^CtJ^p  ) are  long  compared  to  (l/2W),  the  generalized  crosscor- 
relation in  integral  form  can  be  'written  as  follows: 


/:/: 


TW„ 


™ B 1 w, 

a(t-u)v(u)S  (t;cp^)  dudt  = — p 

° ?W  

Q P,q=i 


V a S®  (tp  ) + Vp  S®  (cp  )1 

/ V pqL  IP  W o 2p  iiq  oj 

(2.47) 

The  next  step  is  to  obtain  the  generalized  crosscorrelation 
in  the  form 


TF 


m,  n=l 


a®  V S®(cp  ) 
ran  m n o 


where 


mn 


is  a sample  of  the  kernel  a(t-u)  at  the  times  t^  - m/F 


u = n/F,  and 


I 

I 


V is  a sample  of  V(u)  at  the  time  u = m/F,  etc. 


Ttl 


The  crosscorrelation  L X(t)S^(t;cp^)  dt,  where  X(t)  is  defined  by 
Eq.  (2.30),  is  also  equivalent  to  an  operation  on  the  evenly  spaced 
samples  of  X(t)  and  S®(t;cp  ) obtained  at  a sampling  rate  of  F;  namely. 


m 


£ 


X(t)S®(t;cp^)  dt  = - 


X„sf(cp^) 


(2.48) 


n=-» 


r 

i= 

! 


12 


C.  Van  Schooneveld,  "Some  Remarks  on  Sampling  Methods  for  a Band 
Pass  Signal,"  Signal  Processing,  NATO  Advanced  Study  Institute, 
(September  1964),  p.  413. 


51 


i i 


where 


t = n/F,  n = 0,  ±1,  ±2,  ...  and  X(t  ) = X , etc.  The 
n ' n n 

sampling  rate  F is  sufficiently  high  such  that  X(t)  and  S (t;cp^) 
can  be  reconstructed  by  an  interpolation  formula  from  their 
respective  samples. 

The  sample  X(t  ),  defined  by  Eq.  (2.30),  can  also  be  written 


in  tenns  of  evenly  spaced  samples  of  a(t^  - u)  V(u)  obtained  at  the 


rate  F;  namely, 


X(t  ) = X = i y a^(t  - V(^)  = i y a"  V 

n n F n F F F mn  i 

m=-oo 


mx  ,,/mx 


m 


(2.i^9) 


m=-«> 


where 


a^  = a^  and  V(^)  = V , etc. 
nm  Tim  F m 


Substitution  of  Eq.  (2.49)  into  Eq.  (2.48)  results  in 

/OO  00 


X(t)S  (t;cp^)dt  = A 
F 


O'  „2  2 


V S®(cp  ) 
mn  m n o' 


(2.50) 

niii  Ui  ii 

in,  n=-'» 

Upon  the  application  of  similar  arguments  that  showed  the 
generalized  crosscorrelation  is  approximately  equal  to  am  operation 


on  TW  samples  of  each  quadrature  component,  the  generalized  cross- 


correlation is  also  approximatley  equal  to  an  operation  on  TF  samples 


■g 

each  of  V(u)  and  S (t;cp^);  namely. 


i(t 


u)  V(u)S  (t;cp  )dudt  = — jr  a'^  'V  S^(cp  ) 

o „2  mn  m n o 


TF 


m,  n=  1 


(2.51) 


J 
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substitution  of  Eq.  (2.47)  into  Eq.  (2.5i)  results  in 


V S®(cp  ) = ~^y  a 
2 ^ mn  m n'^o  ^2  P' 

m,n=l  Q p, q=l 


a '^fv  S®  (cp  ) + V S® 
pqL  Ip  Iq'^o^  2p  2q 


<”o>]  ’ 


(2.52) 


where 


a is  a sample  of  the  kernel  a(t-u)  at  the  times  t = ^ , u = and 
mn  ^ m F n F' 

W 

a is  a sample  of  the  kernel  a(t-u)  at  the  times  t = ^ , 
pq  ^ 

u = , etc. 

'1  Wq 

The  function  S (t;cp^)  can  be  written  as 

S®(t;cp  ) = S®(t;cp  ) cos  co  t + S®(t;cp  ) sin  co  t,  or 
o 1 o o 2 o ’ 

S^(t;q5^)  = S®(t+e^)  cos  o:i^(t+e^)  + S®(t+e^)  sin  cn^(t+G_^)  , 


where 


If  -n  < cp^  < jt,  then  it  is  possible  to  make  the  approximations 
S®(t+£o)  = S®(t)  and  S®(t+e^)  s S®(t),  and  S(t;cp^)  becomes 

S®(t;cp  ) = S®(t)  [cos  CD  t cos  cp  - sin  cd  t sin  cp  1 
o 1 L o o o oj 

+ S®(t)fsin  CD  t cos  (p  + cos  cd  t sin  <p  1 
2 L o o o <q) 

If  S®(t;cp^)  is  sampled  at  the  times  t = » q = 0,  ±1,  ±2,  ...,  then 

o q 


f -t) 

®lq^^Po^  " ®lq  'Po  ®fq  'Po 


sin  cp  , or 

o ' 


..  'I*'-.  . 


Mipi  »i  <4 


(2.55) 
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B 

If  S (t;cp^)  is  sampled  at  the  times  t = ^ + , q = 0,  ±L,  ±2, 

then 


_B/^qK  ,1  \ ^s/qK  1 \ ^b(c 

\r  ^ w ’'Poj " ^ irr 


. aK 


f ’^O, 
o 


= -®lq  ^ ^q  'Pq  ' (2.5^^) 

Substitution  of  Eq,  (2.53)  and  Eq,  (2.54)  into  Eq.  (2.52) 


results  in 


TW,, 


W. 


TF  Q „ 

L,  J V s%)  = ^ a ^ is®  cos  cp  + S®  sin  cp 

^ mn  m n 2 / pqL  Ip  \ Iq  ^o  2q  ^o 

m,n=l  Q P,q=4 


+ (-S®  sin  cp  + S®  cos  cp 

2p  \ Iq  2q  o 


, or 


TF 


F „B 


^ a V S 
f -•  mn  m 

m,  n=l 


n^^)=  -y, 

2W„  ^ 


P w 

F"  ^ 

3« 


pq 


Q p,q=I 


V,  S®  + S®  Icos  cp 
Ip  Iq  tip  2q' 


2 ™q  w 

. jL.  v'^  y 


(2.55) 


2W, 


Q P,q=i 


V,  S®  - V.  S®  I sin  cp 
pq|  Ip  2q  2p  Iq' 


Upon  the  application  of  similar  arguments  that  resulted 

in  Eq.  (2.55)j  it  follows  that 

TW 

y af  S®(<p  )s®(cp  ) . ^ Y'‘‘  a“«rs=  S®  a s®  S?  1 , {2.56) 

Zi^  mn  m o n o ^^2  / j pql_  Im  In  2m  ^nj  ' v > / 


m,n=l 


Q p,q=l 


.2 


Y / V S%  ) = I.  y'  a ^[v  + V Icos  cp 

mn  m n o pq[_lp  Iq  2p  2qJ  o 


m,n=l 


^ p,q=l 

TW_ 

p2  Q 


y a 

2w"  P'lL 


V,  1 sin  cp 

Ip  2q  dp  IqJ 


, and 


Q p,q=i 


TF 

E ‘ 

Tti,n=l 


2 '^0  W r 1 

S-^(rp  )S^(cp  ) = ^ y a ^ 
ran  pqL  Lm  In  an  2nJ 


Q p,q=l 


In  order  to  shorten  the  equation  describing  A(V)  in  terms 
of  quadrature  components,  the  following  substitutions  are  made. 


c/  . jf  ■ 

,2  2Li  pq[  Ip  Iq  2p  2q 

Q p,q=l 

TW 

a“‘»rs®  S®  S®  S® 
k„2  pqL^P  2p  ^qj 


Q p,q=t 

TW 

p2  3 ^Qf  ^ M 

= V a ^ V + v„  S 

2^2  Z-  pq[  IP  Iq  2p  2q 


Q p,q=l 


Q TW 

,2  q W 


_I_  « ^fv  S®  + V S®  ' 

“ 2w^  ^ pq|_  ip  iq  2p  2q_ 


Q.  p,q=i 


— V - V 

" 2 ^ pq  L ip  2q  ■ 2p  Iq 


Q p,q=i 


TW 

2 0 W n 

2 2^  pq  Ip  2q  2p  Iq 

0.  P,q=i 


6m  • 


(2.64) 
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The  equation  for  A(V)  then  becomes 


A(V)  = e 


/: 

/ 


B B 

« \ cos  cp  +v  sin  CD 

A A 

■t  \ cos  cp  +v  sin  cp 

O C5  , 

e dcp 


(2.65) 


If  the  following  change  of  variables  is  made,  'T’q  " 


the  result  is 


A(^) 


I: 

I 


r-  B B 

'~2n  cos  cp'  - V sin  cp' 

o ° 

cZ-c^  Jo  ' 


A A 

2n  COE  cp'  - V sin  cp' 


dcp' 

o 


The  integrals  can  be  integrated  directly  to  give 
A(^)  = 


I 

0 

7 _2  _2^ 

K -v®j 

|i/3 

7 /,2 

f,  A A 

2- 

I 

0 

+ V J 

(2.66) 


(2.67) 


where 

I represents  the  modified  Bessel  fionction  of  the  first  kind  of  order 
o 


zero. 


The  equation  for  i(V),  from  Eq.  (2.22),  is 

- y'  a^  (v  - V S®  i S®S®1 

/ .■  mnlmn  mn2mn  2mnJ 

£(V)  = e .(2.68) 

Written  in  terms  of  quadrature  components,  Eq.  (2.68)  becc^mes 


-2 


2W^  P<lL 

= e ^ 


A A R B 

V,  S,  + V„  - V, 

Ip  Iq  2p  2q  Ip  Iq  2p  2q 


(2.69) 


_!  qA  A i c-^  + i S®  S®  + i S®  S® 

■ 2 ®lp®lq‘  2 ®2p®2q  2 ®lp®lq  2 2p®2q_ 
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One  advantage  in  using  quadrature  components  for  likelihood 


ratio  calculations  is  that  the  dimensions  of  the  square  matrix 


are  smaller  than  the  dimensions  of  the  square  matrix 


faM. 

\ mn/ 


(.y 


A second 


advantage  is  that  the  statistics  of  the  in-phase  components  fmd  the 


out-of-phase  components  of  Gaussian  noise  remain  Gaussian. 


CHAPTER  III 


DATA  USED  FOR  LIKELIHOOD  RATIO  PROCESSING 


A.  Arrangement  of  Apparatus 

The  data  were  taken  at  the  Lake  Travis  Test  Facility  of  the 
Applied  Research  Laboratories  of  The  University  of  Texas  at  Austin, 
The  source  transducer,  target,  and  receiving  hydrophone  were 
suspended  in  a line  beneath  a barge  at  a depth  of  l6  feet,  as  in- 
dicated in  Fig.  3.1.  The  source  transducer  (Edo  327)  was  in  effect 
a "piston"  transmitter  with  a resonant  frequency  of  70  kHz . 

The  output  Y(t)  of  a transducer,  given  the  input  X(t)  and 
the  impulse  response  H(t)  of  the  transducer,  is 


Y(t)  = / X(t-T)  H(x)dT 


(5.1) 


For  an  input  signal  bandlimited  to  the  frequency  interval 
f^  - W/2  < |f|  < f ^ + w/2,  Y(t)  can  be  written  as 


r-f  4 

■/,v 


X(f)H(f)  e’-^^^df  + 


f .Vi 

f 


X(f)H(f)  e^^^df  . (3.2) 


H(f)  in  general  is  a complex  function,  and  it  can  be 


written  as 


H(f)  = A(f)  , (3.3) 

where  A(f)  is  the  amplitude  distortion  and  cp(f)  is  the  phase 
distortion  of  the  transducer. 

For  the  system  composed  of  the  source  and  receiving  trans- 
ducers, A(f)  was  "flat"  in  the  frequency  interval  of  60  kHz  to  80  kHz, 


■4 


i RECEIVING  HYDROPHONE 


V 

SOURCE  TRANSDUCER 


TARGET 


FIGURE  3.1 

DIAGRAM  OF  THE  EXPERIMENTAL  ARRANGEMENT 
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even  though  the  source  transducer  had  a resonant  frequency  at  JO  kHz . 
(The  amplitude  distortion  is  considered  "flat"  in  a frequency  interval 
if  there  is  less  than  a 5 dB  change  in  the  magnitude  of  the  amplitude 
distortion  in  the  frequency  interval. ) Since  the  source  transducer 
had  a resonant  frequency  at  JO  kHz,  optimum  use  was  made  of  this 
sensitivity.  For  the  system,  cp(f)  was  linear  in  the  frequency  interval 
of  60  kHz  to  80  kHz.  (The  linearity  of  cp(f)  in  a frequency  interval 
results  only  in  a time  delay  of  a signal  in  the  time  domain.) 

The  targets  were  suspended  from  the  barge  by  0.015  in.  diam 
stainless  steel  wire,  22  feet  from  the  source  transducer.  The  re- 
ceiving hydrophone  was  located  between  the  source  and  the  target,  l6 
feet  from  the  target. 

B.  Targets 

The  five  targets  used  were 

(1)  a 7 in.  diam  solid  aluminum  sphere, 

(2)  a 5 in.  diam  solid  aluminum  sphere, 

(5)  a 5 in.  diam  solid  aluminum  sphere, 

(4-)  a 5 in.  diam  styrofoam  sphere,  and 

(5)  a 5 in.  diam  hollow  aluminum  sphere  with  a b/a  of  0.9T5. 

(b  is  the  inner  diameter  of  the  hollow  sphere  and  a is  the  outer 
diameter.)  The  region  inside  the  hollow  sphere  contained  air  at  a 
pressure  of  one  atmosphere. 

The  two  signal  transmissions  used  were 

(l)  a sinusoidal  pulse  of  frequency  JO  kHz  with  a pulse 


length  of  80  usee,  and 


ko 

(2)  a linear  M of  bandwidth  12.5  kHz  with  a center 
frequency  of  65.75  Wiz  and  a pulse  length  of  3 msec. 

Plots  of  the  two  signal  transmissions  are  shown  in  Fig.  3.2.  In 
Table  3.1,  ka  values  for  a 3 in.  diam,  a 5 in.  diam,  and  a 7 in.  diam 
sphere  for  the  2 different  transmissions  are  listed. 

C . Echoes  from  the  Targets 

"The  echoes  from  a rigid,  freely  movable  sphere  differ 

significantly  from  those  from  an  immovable  sphere  only  for  values  of 

ka  less  than  flve."^^  For  values  of  ka  less  than  two,  the  echoes  from 

the  rigid,  freely  movable  sphere  are  similar  to  those  from  the  solid 

elastic  sphere  for  corresponding  values  of  the  densities  of  the  two 

spheres. "Presumably,  beyond  this  stage,  elastic  vibrations  have  an 

increasing  effect  until  they  become  the  dominant  factor  in  determining 

the  very  marked  frequency  dependence  of  echoes  from  solid,  elastic 

15 

materials  in  water." 

For  the  ka  range  of  10.5  to  26.9,  an  echo  from  a completely 
rigid  sphere  is  similar  to  that  of  the  incident  pulse,  and  the  arrival 
time  of  the  echo  shows  that  it  was  reflected  from  the  portion  of  the 
the  sphere  closest  to  the  source. The  echoes  from  elastic  spheres, 
such  as  the  solid  aluminum  spheres,  are  longer  and  more  ccmplicated . 


13 

Robert  Hickling  and  N.  M.  Wang,  "Scattering  of  Sound  by  a Rigid 
Movable  Sphere,"  J.  Acoust.  Soc.  Am,  Ho,  2 (1966)  p.  276. 

^^Ibid.,  p.  276. 

^^Ibid.,  p.  276. 

^ K.  Jerome  Dierks  and  R.  Hickling,  "Echoes  from  Hollow  Aluminum 
Spheres  in  Water,"  J.  Acoust.  Soc.  Am.  No.  2 (1967)  p.  38O, 
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LINEAR  FM  TRANSMISSION 

f = 63.75  kHz  W = 12.5  kHz 
o 

T = 3 msec  At  = 1 ^300  msec 


: 70  kHz  T = 80  fisec 
At  = 1 300  msec 

i 

I FIGURE  3.2 

’ TRANSMISSION  WAVEFORMS 

PEAK  |AMPLITUDE|  = 0.5  in. 
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The  first  portions  of  the  echoes  still  retain  the  characteristics  of 

the  rigid  sphere  echoes,  indicating  that  there  are  primary  reflections 

from  the  outer  surface,  similar  to  those  from  a rigid  body,  followed 

17 

by  additional  echoes  caused  by  the  elastic  reaction  of  the  sphere. 

When  the  length  of  the  incident  pulse  is  less  than  the  diameter  of 

the  sphere,  the  primary  echo  is  separated  in  time  from  the  secondary 
l8 

elastic  echoes.  As  the  incident  pulse  length  becomes  shorter  (or  as 

the  diameter  of  the  sphere  increases),  each  of  the  component  parts  of 

19 

the  echo  contract  and  separate.  As  the  incident  pulse  grows  longer 

(or  as  the  diameter  of  the  sphere  decreases),  the  secondary  echoes, 

due  to  the  elastic  vibrations  in  the  sphere,  begin  to  overlap  with  the 

primary  reflection  and  with  each  other,  and  the  interactions  interfere 

constructively  or  destructively,  depending  on  the  particular  value  of 

20 

the  dominant  frequency. 

The  first  highlights  of  the  echoes  from  a 5 in.  diam  hollow 

aluminum  sphere  with  b/a  = 0.975  correspond  to  rigid  body  reflections. 

The  remaining  portions  of  the  echoes  are  a result  of  the  vibrations 

21 

induced  within  the  sphere.  For  short  incident  pulses,  these 
secondary  echoes  are  separated  from  the  primary  echoes . For  long 
pulses,  the  secondary  echoes  overlap  and  interfere  with  the 


^'"ibld.,  p.  581. 

^^^Ibid.,  p,  381. 

^^Ibid.,  p.  381. 

^^Ibld.,  p.  381, 

21 

Robert  Hickling,  "Analysis  of  Echoes  from  a Hollow  Metallic  Sphere 
in  Water,"  J.  Acoust,  Soc . Am.  No.  6 (1964)  p.  1153. 
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primary  echo.  For  b/a  = 0.975,  the  waves  that  occur  in  the  solid 

22 

material  of  the  shell  are  mainly  of  a flexural  type.  These  waves 

travel  around  the  shell  with  a velocity  the  same  as  or  slightly  less 

25 

than  the  velocity  of  shear  waves  in  the  material. 

The  5 in.  diam  styrofoam  sphere  can  be  approximated  as  a 

5 in.  diam  sphere  with  a pressure  release  surface.  For  a ka  of  about 

20,  the  reflected  pulse  has  the  same  wave  shape  as  the  emitted  pulse, 

24 

but  the  reflected  pulse  is  inverted. 

The  echoes  for  the  sinusoidal  pulse  transmission  were 
passed  through  a bandpass  filter  of  bandwidth  (w)  equal  to  20  kHz 
and  a center  frequency  (f^)  equal  to  70  kHz.  The  filter  was  designed 
to  eliminate  frequencies  outside  the  20  kHz  bandwidth  without  the 
distortion  of  frequencies  within  the  12.5  kHz  bandwidth.  The  echoes 
for  the  linear  FM  transmission  were  passed  through  a bandpass  filter 
of  bandwidth  (w)  equal  to  20  kHz  with  a center  frequency  (f^)  equal 
to  65  kHz . Again  the  filter  was  designed  to  eliminate  frequencies 
outside  the  20  kHz  bandwidth  without  the  distortion  of  frequencies 
within  the  12.5  kHz  bandwidth.  The  following  data  were  recorded  on 
magnetic  tape.  On  one  channel  of  the  tape,  the  transmitted  signal 
was  recorded.  On  a second  channel,  a 500  kHz  continuous  reference 


^^Ibid.,  p.  1155. 
^^Ibid . , p.  1155. 


Robert  Hickling,  "Analysis  of  Echoes  from  a Solid  Elastic  Sphere 
in  Water,"  J.  Acoust.  Soc.  Am.  No.  10  (1962)  p.  1584. 


tone  was  recorded,  and  on  the  third  channel,  the  echoes  from  the 
spheres  were  recorded.  The  300  kHz  reference  tone  was  used  in  the 
digitizing  process  to  eliminate  the  "wow  and  flutter"  of  the  tape 
recorder.  The  transmitted  signal  was  used  as  a trigger  for  the 
analog-to-digital  converter  to  begin  the  analog-to-digital 
conversion  of  an  echo. 

The  echoes  were  digitized  at  a rate  of  300  kHz  and  recorded 
in  digital  form  on  magnetic  tape.  The  gain  settings  of  the  amplifiers 
used  in  the  digitizing  process  were  set  so  that  the  amplitudes  of  the 
analog  echoes  had  a possible  maximum  absolute  value  of  10  volts.  In 
the  analog-to-digital  conversion,  the  amplitudes  of  the  digitized 
echoes  were  scaled  so  that  10  volts  were  equal  to  20^7 . The  number 
20kT , (2^^  -l),  resulted  from  digitizing  the  amplitudes  of  the  echoes 
into  12  bit  words  with  the  first  bit  of  each  word  used  as  a sign  bit 
(12  bit  quantization). 

In  Fig,  3.3  an  echo  from  each  target  for  a sinusoidal  pulse 
transmission  is  shown.  The  leading  edge  of  the  envelope  for  each 
echo  in  Fig,  3.3  is  characterized  by  a steep,  positive  slope. 

In  Fig,  3.^  an  echo  for  a linear  M transmission  from  each 
target  is  shown.  The  echoes  for  a linear  FM  transmission  were  passed 
through  a "gate"  for  a duration  of  5 msec.  The  gate  was  used  to 
minimize  the  surface  and  volume  reverberation  in  the  recorded 


returns . 
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FIGURE  3.3 

; > . SONAR  RETURNS  FOR  A SINUSOIDAL  PULSE  TRANSMISSION 

PULSE  LENGTH  s 0.37  msec  EPOCH  = At  = 1/300  msec 

f = 70  kHz  W = 20kHz  PEAK  IaMPLITUDEI  : 0.5  in. 
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ECHO  FROM  5 in.  diam  HOLLOW  ALUMINUM  SPHERE 


ECHO  FROM  5 in.  diom  STYROFOAM  SPHERE 


! ECHO  FROM  3 in.  diam  SOLID  ALUMINUM  SPHERE 


ECHO  FROM  5 in.  diam  SOLID  ALUMINUM  SPHERE 


ECHO  FROM  7 in.  diam  SOLID  ALUMINUM  SPHERE 


FIGURE  3.4 

SONAR  RETURNS  FOR  A LINEAR  FM  TRANSMISSION 
PULSE  LENGTH  = 3.2  msec  EPOCH  - At  - 1 '300  msec 
s 65  kHz  W=  20  kHz  PEAK  |AMPLITUDE|  = 0.5  in. 
DURATION  OF  GATE  = 5 msec 


CHAPTER  IV 


TARGET  CLASSIFICATION  USING  A 
SIIWSOIDAL  PULSE  TRANSMISSION 


A.  Classification  of  Simllax  Targets 

1.  Signal-to-Noise  Ratio 

The  ability  of  the  likelihood  ratio  processor  to  classify 
sonar  returns  from  similar  targets  for  a sinusoidal  pulse  transmission 
was  investigated.  The  five  targets  used  were  the  J in.  dlam, 

5 in.  diam,  and  3 in.  diam  solid  aluminum  spheres,  the  5 in.  diam 
styrofoam  sphere,  and  the  5 in.  diam  hollow  aluminim  sphere.  The 
sinusoidal  pulse  transmission  and  the  echoes  from  the  targets  were 
described  in  Chapter  III. 

The  s/N  (signal-to-noise)  ratio  of  an  echo  was  defined  as 

geak  _ynplltude  of  signal  ^ units  of  decibels  (dB),  the  S/N  ratio 
rms  of  noise  amplitude  \ /y  / 

, / peak  amplitude  of  signal\ 

was  expressed  as  20  log^_  ^ . rrr— § dB- 

10  '\^rms  of  noise  amplitude  / 

The  S/N  ratios  of  the  digitized  echoes  from  the  five  targets 
were  large,  ranging  from  an  average  S/n  ratio  of  3^  dB  for  the  echoes 
from  the  3 diam  solid  aluminum  sphere  to  an  average  S/N  ratio  of 
4l  dB  for  the  echoes  from  the  5 in.  diam  hollow  alumimma  sphere. 

2 . Epoch 

The  use  of  a threshold  technique  to  detemnine  the  epoch 


G = — j of'  4he  digitized  echoes  from  the  five  targets  was  possible 
o 0)^ 

since 


(1)  the  s/n  ratios  of  the  echoes  were  large, 

(2)  the  variations  in  A(t)  and  T(t)  from  echo  to  echo  among 
the  echoes  from  the  same  target  were  small,  and 
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1+9 

(3)  the  leading  edge  of  the  envelope  for  each  echo  was 
characterized  by  a steep  positive  slope. 

This  steepness  in  the  slope  was  an  important  factor  in  the  determina- 
tion of  epoch.  The  epoch  (e^)  for  an  echo  was  defined  as  the  time  at 
which  the  amplitude  of  an  echo  reached  20^  of  its  peak  amplitude. 

After  epoch  w-as  established,  a new  epoch  (^q)  vas 

defined  in  order  to  include  in  a likelihood  ratio  calculation  the 
portion  of  an  echo  occurring  before  the  threshold.  The  new  epoch  (e^) 
was  defined  as  occurring  0.027  msec  (or  8 samples)  before  the  threshold. 
The  pulse  lengths  of  the  echoes  were  defined  to  be  O.JT  msec,  as 
indicated  in  Fig. 

In  Fig.  l+.l,  the  success  of  a threshold  technique  in 
determining  epoch  (e^)  is  demonstrated.  Seven  echoes  from  a 5 in.  diam 
solid  aluminum  sphere  are  shown.  Each  echo  is  plotted  along  a different 
time  axis.  The  distance  from  the  beginning  of  a time  axis  to  the  epoch 
is  the  same  for  each  echo. 

The  amplitudes  of  each  echo  from  each  target  were  scaled  so 

that  the  peak  amplitude  of  each  echo  from  each  target  was  equal  to  l400. 

The  purpose  of  this  scaling  was  to  be  able  to  base  the  classification 

of  targets  on  the  differences  in  waveform  of  the  echoes.  The  number 

of  echoes  available  from  each  target  varied  from  12  to  l4.  A Hffer- 

ence  echo  for  a target  was  computed  by  taking  an  ensemble  average  of 

all  the  echoes  from  the  target.  The  reference  echo  for  target  "A  was 

defined  as  > , and  the  reference  echo  for  target  B was 

K 

defined  as  > . (The  brackets  denote  an  ensemble  average 

K 

over  the  available  set  of  target  echoes.) 
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£q|« — 0 37  mMc  -n 


FIGURE  4.1 

ECHOES  FROM  A 5 in.  diam  SOLID  ALUMINUM  SPHERE 

PULSE  LENGTH  = 0.37  msec  EPOCH  = 

At  = 1 /300  msec  = 70  kHz  W --  20  kHz 
= 38.8  dB  PEAK  |AMPLITUDE|  = 0.5  in. 


3.  Addition  of  Gaussian  Noise  to  Echoes 

Samples  of  stationary  Gaussian  noise  of  mean  zero  from  a 
random  noise  generator  were  added  to  all  the  echoes  for  two  purposes : 
(l)  to  make  the  simulation  of  target  classification  more 


realistic; 

(2)  to  increase  the  number  of  messages  (echo-noise 
combinations)  from  each  target  class.  (The  members  of  a target  class 
are  characterized  by  the  same  Gaussian  noise  statistics  and  the  same 
reference  echo.) 

The  bandwidth  of  the  Gaussian  noise  was  20  kHz  with  a center  frequency 

of  70  kHz.  A power  spectrum  computed  from  samples  of  the  Gaussian  noise 

25 

is  plotted  in  Fig.  4.2.  The  confidence  interval  used  in  pl...o-l;;£  the 
power  spectrum  was  8o5.  The  random  noise  generator  is  further  discussed 
in  Appeiidix  A. 

The  problem  of  classification  was  made  more  realistic  by 
decreasing  the  s/N  ratios  of  the  echoes  from  a target  to  a prescribed 
lower  value  uj  si^ecifying  the  rms  of  the  amplitude  of  the  Gaussian 
noise  to  be  added  to  the  echoes  from  the  target.  The  desired  S/N  ratio 
of  the  echoes  from  a target  was  described  as  follows: 


(4.1) 


< > + < N^  > 


'^R.  B.  Blackman,  and  J.  W-  Tukey,  The  Measurement  of  Power  Spectra.. 
Dover  Publications,  Inc.,  New  York,  N.  Y. , (1958),  p.  23. 
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where 
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= lUOO  = peak  amplitude  of  each  echo  from  a target, 


<4 

< 


> = mean  square  of  ambient  noise  pi-esent  in  the  echoes, 

> = mean  square  of  noise  from  the  random  noise 

generator,  and 

= desired  S/N  ratio  of  each  echo  from  a target. 


The  ambient  noise  in  an  echo  was  independent  of  the  additive  Gaussian 
noise  from  the  random  noise  generator.  For  of  the  order  1,  the 
rms  of  the  amplitude  of  the  ambient  noise  in  an  echo  was  much  less  than 
the  rms  of  the  amplitude  of  the  Gaussian  noise  from  the  random  noise 
generator.  From  Eq.  (4.1)  the  rms  of  the  amplitude  of  the  noise  from 
the  random  noise  generator  was  specified  as  follows: 


(‘•.2) 


For  Rgjj  of  the  order  1,  the  rms  of  the  amplitude  of  the  ambient  noise 

was  neglected  in  the  computation  of  < > 

In  the  investigation  of  the  ability  of  the  likelihood  ratio 
processor  to  classify  similar  targets,  Rg^j  was  set  equal  to  one  (O  dB). 
In  Fig.  4.5,  a message  with  a S/n  ratio  of  0 dB  from  each  target  class 
is  shown.  In  Fig.  4. 5,  the  distorted  echoes  of  pulse  length  0.37  msec 
are  located  between  the  two  vertical  lines. 

The  number  of  messages  from  each  target  class  was  increased 
by  the  addition  of  10  separate  noise  signals  from  the  random  noise 
generator  to  each  echo  so  that  the  total  number  of  messages  from  each 
target  class  varied  from  120  to  l40.  No  noise  signal  was  duplicated 
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H*  0 37  msec  “i 


ECHO  FROM  5 in.  diam  HOLLOW  ALUMINUM  SPHERE  PLUS  NOISE 


i 1 

ECHO  FROM  5 in.  diam  STYROFOAM  SPHERE  PLUS  NOISE 


ECHO  FROM  3 in.  diam  SOLID  ALUMINUM  SPHERE  PLUS  NOISE 


» 

ECHO  FROM  5 in.  diam  SOLID  ALUMINUM  SPHERE  PLUS  NOISE 


ECHO  FROM  7 in.  diam  SOLID  ALUMINUM  SPHERE  PLUS  NOISE 


FIGURE  4 3 

MESSAGES  FOR  SINUSOIDAL  PULSE 

POSITION  OF  0.37  msec  TIME  INTERVAL  INDICATES  LOCATION  OF  EMBEDDED  ECHO 
PULSE  LENGTH  = 0.37  msec  At  = 1 '300  msec  = 70  kHz 

W = 20  kHz  S/'N  = 0 dB  PEAK  |AMPLITUDE|  = 0.5  in. 
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throughout  this  procedure.  In  the  classification  of  the  messages 
from  two  target  classes,  the  same  number  of  messages  from  each  target 
class  was  used.  Therefore,  the  a priori  probability  (q^^)  for  the 
messages  from  target  class  A was  equal  to  the  a priori  probability  (q^) 
for  the  messages  from  target  class  B. 

The  knowledge  of  epoch  was  retained  during  the  addition  of 
Gaussian  noise  to  the  echoes  so  that  the  likelihood  ratio  processor 
described  by  Eq.  (2.69)  could  be  used  to  classify  the  messages. 

In  obtaining  the  quadratirre  components  of  the  messages  by 
quadrature  sampling,  a sampling  parameter  of  K = 2 was  used.  This 
sampling  parameter  gave  the  minimum  number  of  samples  (TW')  greater  than 


TW'  of  the  in-phase  component 


= 30  kHz,  f 
’ o 


75  kHz,  and  W’ 


and  the  minimum  number  of  samples  (TW')  greater  than  TW'  of  the 
out-of-phase  component.  This  sampling  parameter  was  used  in  the  com- 
putation of  the  envelope  plotted  in  Fig.  2.5* 

U.  Computation  of  Autocovariance  Matrix  for  Noise 


An  autocovariance  matrix  C of  order  TW'  x TW'  was  computed 
from  a discrete  noise  signal  using  the  following  formula: 

P-T 


N N , 

a a+T  ’ 


(^O) 


(7=1 


where 


26 

P g 10  TW ' , and 

y 


N = N 
a 


(k) 


Julius  S.  Bendat,  and  A.  G.  Piersol,  Measurement  and  Analysis  of 
Random  Data,  John  Wiley  and  Sons,  Inc.,  New  York,  N.  Y. , (19(ib),  p.  29O. 
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The  noise  signal  consisted  of  samples  of  stationary  Gaussian  noise  from 
the  random  noise  generator.  The  rms  of  the  noise  signal  amplitude  was 


< 


1/2 


> 


SN 


(1.4) 


where 


= desired  S/N  ratio  of  each  message,  and 
= l400  = peak  amplitude  of  each  echo. 

The  bandwidth  of  the  noise  signal  was  20  kHz  with  a center  frequency  of 

W'  W'  W ' 

70  kHz.  The  matrix  a ^ of  order  TW^^  x TW^  where  a ^ C = 1,  was  com- 

piited  by  taking  the  inverse  of  the  autocovariance  matrix  C . 

5.  Geometrical  Description  of  Target  Classes 

A measure  of  the  ability  of  the  likelihood  ratio  processor 
to  classify  messages  for  a fixed  threshold  can  be  determined  geomet- 
rically as  follows;  The  message  vectors  from  target  class  A are 
labeled  as  and  the  message  vectors  from  target  class  B are  labeled 
as  The  vectors  A and  are  located  in  a 2TW' -dimensional 

vector  space  and  are  depicted  as  located  in  a 3-<ii™enslonal  vector  space 
in  Fig.  4.4. 

The  geometrical  representation  in  Fig.  4.4  is  based  on  the 
assumption  that  (a  ) = (a  & ) where  a is  a constant.  In  Fig.  4.4, 

the  average  radius  of  the  sphere  generated  by  the  vectors'^  is 

< > , and  the  average  radius  of  the  sphere  generated  by  the 

vectors  is  <\i^  - tl\  > . The  deviation  of  ["v  - from 

K A 

< ""^1  > is  < ^<  1"^  - I > - > , and  the  deviation 


c. 


«. 
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of  -'^®|  from  < It®  I > is  < (<  - -g®  I > - If®  - . 

The  distance  between  the  centers  of  the  two  spheres  is  - t®|  . Tiie 

shaded  portion  indicates  the  overlapping  volume  of  the  two  spheres. 

For  (a  ) (a  6 ),  the  spheres  are  distorted  into  ellipsoids, 

mn  ' o mn  ^ ^ 

A hyperplane  was  constructed  to  intersect  the  centers  of  the 
two  hyperspheres.  A schematic  projection  of  the  hyperspheres  onto  the 
hyperplane  is  depicted  in  Fig.  4.5. 

In  Fig.  4.5,  the  short  lines  perpendicular  to  the  perimeters 
of  the  two  circles  indicate  the  deviation  of  the  radii  of  each  circle 
from  the  average  radius.  If  | r^  - r^  | and  |d^  - d^|  are  small  com- 
pared to  r^  and  d^  respectively,  the  geometric  threshold  for  equal 
costs  for  wrong  decisions  and  equal  a priori  probabilities  would  be  a 
line  that  perpendicularly  bisects  - f®( . On  consideration  of  the 
restrictions  in  the  last  statement,  a geometric  measure  of  the  ability 

®T 

to  classify  two  targets  is  — , where 

o 

D < iv* -th  > + < It® l>+ -^ol 

It  is  possible  for  to  go  to  zero  without  the  attainment 
of  perfect  classification.  An  error  in  classification  can  still  occur 
since  the  radii  of  each  sphere  deviate  from  an  average  radius.  The 
angle  between  the  two  reference  vectors  and  f®  would  be  an  unsatis- 
factory measure  of  the  ability  of  the  likelihood  ratio  processor  to 
classify  targets,  because  it  does  not  indicate  the  values  |"^|  and 


FIGURE  4.5 

INTERSECTION  OF  THE  CENTERS  OF  THE  HYPERSPHERES 
WITH  A HYPERPLANE  REPRESENTED 
IN  A 2-DIMENSIONAL  SPACE 
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6.  ROC  Curves 

Another  measure  of  the  ability  of  the  likelihood  ratio 
processor  to  classify  two  targets  is  an  ROC  curve.  (ROC  is  an  acronym 
for  Receiver  Operator  Characteristic.)  An  ROC  curve  gives  the  proba- 
bility of  correctly  classifying  target  A and  the  probability  of 
incorrectly  classifying  target  B for  any  threshold. 

An  ROC  cu2Tve  is  generated  as  follows:  A threshold  is 

specified.  The  number  of  likelihood  ratios  computed  for  the  messages 
from  target  class  A that  are  less  than  the  threshold  are  counted. 

This  number  is  normalized  by  dividing  it  by  the  total  number  of  like- 
lihood ratios  computed  for  the  messages  from  target  class  A.  The 
resulting  nimiber  is  the  probability  of  correctly  classifying  target  A 
for  the  specified  threshold.  Next,  the  number  of  likelihood  ratios 
computed  for  the  messages  from  target  class  B that  are  less  than  the 
threshold  are  counted.  This  number  is  normalized  by  dividing  it  by 
the  total  number  of  likelihood  ratios  computed  for  the  messages  from 
target  class  B.  The  resulting  number  is  the  probability  of  incorrectly 
classifying  target  B for  the  specified  threshold.  The  above  procedure 
is  repeated  until  all  possible  values  of  the  threshold  are  used.  For 

‘^1*"21 

zero  costs  for  right  decisions,  the  optimum  threshold  is  — = — , as 

discussed  in  Chapter  II,  Section  A.  The  area  under  an  ROC  curve  gives 

a measure  of  the  performance  of  the  likelihood  ratio  processor  integrated 

over  all  values  of  the  threshold.  The  area  times  100  is  equal  to  the 

27 

expected  percentage  of  correct  decisions. 


27. 


David  M.  Green,  and  John  A.  Swets,  Signal  Detection  and  Psychophysics, 
John  Wiley  and  Sons,  Inc.,  New  York,  N.  Y. , (1966),  p.  217. 
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The  ROC  curves  will  be  plotted  within  a square  frame  and  the 
diagonal  that  bisects  the  lower  left  corner  of  the  frame  represents  a 
performance  equal  to  chance. 

The  ROC  curves  in  Figs.  4.6  through  4.10  were  computed  using 
the  techniques  described  in  this  chapter  and  the  test  function 


TW, 


Q 


a V S®  + V S®  - V - V , 

pq'^  Ip  iq  2p  2q  Ip  Iq  2p  2qJ 


m,n=l 

The  calculation  of  the  test  function  instead  of  the  likelihood  ratio 
described  in  Eq.  (2.69)  does  not  change  the  optimxm  character  of  an 
ROC  cixrve  since  the  test  function  is  a monotonically  increasing  function 
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of  the  likelihood  ratio.  The  areas  under  the  ROC  curves,  along  with 


'■•a  - '■b' 


*A-  V 


, and  — , are  listed  in  Tables  4.1  through  4.5. 
o 


A A 

7.  Geometrical  Interpretation  of  Results 

The  quantities  r^  and  r^  are  dependent  on  the  number  of 
samples  used  to  describe  an  echo  and  on  the  magnitude  of  the  scaling  of 
amplitudes  of  the  echoes. 

A change  in  r^  with  each  target  B is  a result  of  the  random 
nature  of  the  amplitude  of  the  noise  from  the  random  noise  generator. 

A change  in  r with  each  target  B is  a result  of  the  random  nature  of 
the  amplitude  of  the  noise  from  the  random  noise  generator,  plus  the 
use  of  a different  target. 


2^ 


David  Middleton,  Topics  in  Coimunication  Theory,  McGraw-Hill  Book 
Company,  New  York,  N.  Y. , (I965),  p.  25- 
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PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  TARGET  B 

5 in.  diam  SOLID  ALUMINUM  SPHERE 
3 in.  diam  SOLID  ALUMINUM  SPHERE 
5 in.  diam  STYROFOAM  SPHERE 
5 in.  diam  HOLLOW  ALUMINUM  SPHERE 


FIGURE  4.6 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  SIMILAR  TARGETS 
TARGET  A * 7 in.  diam  SOLID  ALUMINUM  SPHERE 


i » 


TABLE  4.1 


TARGET  A = T IN.  DIAM  SOLID  ALUMINUM  SPHERE,  S/N  = 0 dB 


5 in.  diam  3 in.  diam  5 in.  diam  5 in.  diam 

Target  B Solid  Aluminum  Solid  Aluminum  Styrofoam  Hollow  Aluminum 
Sphere  Sphere  Sphere  Sphere 


"a-^bI 


0.009 


0.001  0.001 


O.OOU 


I 

A B 


Area  under 


ROC  Curve 


0.760 


0.782  0.816 


0.789 


PROBABILITY  OF  CORRECTLY  CLASSIFYING  TARGET 


PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  TARGET  B 

7 in  diam  SOLID  ALUMINUM  SPHERE 

3 in  diam  SOLID  ALUMINUM  SPHERE 

5 in.  diam  STYROFOAM  SPHERE 

5 in.  diam  HOLLOW  ALUMINUM  SPHERE 


FIGURE  4.7 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  SIMILAR  TARGETS 
TARGET  A * 5 in.  diam  SOLID  ALUMINUM  SPHERE 


TABLE  4.2 


TARGET  A = 5 IN.  DIAM  SOLID  ALUMINUM  SPHERE,  S/N 


Target  B 


r -r 
A B 


'^A-Sl 


Area  tinder 
ROC  Curve 


7 in . diam 
Solid  Aluminum 
Sphere 

5 in . diam 
Solid  Aluminum 
Sphere 

5 in.  diam 
Styrofoam 
Sphere 

5 in . diam 
Hollov*  Aluminum 
Sphere 

0.009 

c.oou 

0.004 

0.003 

0.06 

0.02 

0.04 

0.05 

1.29 

1.25 

1.30 

1.32 

0.760 

0.719 

j.763 

0.773 

S/N  = 0 dB 


0.2  0.4  0.6  0.8 

PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  TARGET  B 

7 m.  diam  SOLID  ALUMINUM  SPHERE 

5 in.  diam  SOLID  ALUMINUM  SPHERE 

5 in.  diam  STYROFOAM  SPHERE 

5 in.  diam  HOLLOW  ALUMINUM  SPHERE 


FIGURE  4.8 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  SIMILAR  TARGETS 
TARGET  A * 3 in.  diam  SOLID  ALUMINUM  SPHERE 


TABLE  4.3 


TARGET  A = 3 IN,  DIAM  SOLID  ALUMINUM  SPHERE,  S/N 


Target  B 


A 3| 


r. 


VSi 


Area  under 
ROC  Curve 


7 in . diara 
Solid  Aluminm 
Sphere 

5 in . diam 
Solid  Aluminum 
Sphere 

5 in . diam 
Styrofoam 
Sphere 

5 in.  diam 
Hollow  Aluminum 
Sphere 

0.001 

0.004 

0.003 

0.002 

0.05 

0.02 

0.04 

0.03 

1.50 

1.25 

1.24 

1.25 

0.782 

0.719 

0.706 

0.706 

S/N  = 0 dB 


0.2  0.4  0.6  0.8 

PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  TARGET  B 


7 in.  diam 
5 in.  diam 
3 in.  diam 
5 in.  diam 


SOLID  ALUMINUM  SPHERE 
SOLID  ALUMINUM  SPHERE 
SOLID  ALUMINUM  SPHERE 
HOLLOW  ALUMINUM  SPHERE 


FIGURE  4.9 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  SIMILAR  TARGETS 

TARGET  A * 5 in.  diam  STYROFOAM  SPHERE 
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7 in.  diam  5 in.  diatn  5 in.  diam  5 in.  diatn 

Target  B Solid  Alumin\im  Solid  Aluminum  Solid  Aluminum  Hollow  Aluminum 
Sphere  Sphere  Sphere  Sphere 


IV^bI 


0.001 


0.004 


0.005 


IV^I 


0.004 


Area  under 


ROC  Curve 


0.8l6 


0.763 


0.706 


0.581 


I 

I 
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CODE  FOR  TARGET  B 

7 in.  diom  SOLID  ALUMINUM  SPHERE 

5 in.  diom  SOLID  ALUMINUM  SPHERE 

3 in.  dkam  SOLID  ALUMINUM  SPHERE 

5 in.  diom  STYROFOAM  SPHERE 


FIGURE  4.10 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  SIMILAR  TARGETS 
TARGET  A » 5 in.  diam  HOLLOW  ALUMINUM  SPHERE 


i f 
i ' 
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In  Tables  4.1  through  4,5,  (r^  - r^j  and  |d^  - dg(  are  small 
compared  to  r^  and  d^  respectively.  Therefore,  a threshold  of  1 can 
be  interpreted  geometrically  as  a hyperplane  which  perpendicularly 
bisects  the  line  that  connects  the  centers  of  the  two  overlapping 
hyperspheres  generated  by  the  vectors  ^ and 

The  test  functions  computed  from  the  vectors  and  the  test 
functions  computed  from  the  vectors  both  have  Gaussian  distributions. 
The  Gaussian  distributions  determine  the  shape  of  the  ROC  curves  in 
Figs.  4.6  through  4.10. 

All  the  computations  described  in  this  chapter  and  in  later 
chapters  were  done  on  the  CDC  3200  computer  at  the  Applied  Research 
Laboratories  of  The  University  of  Texas  at  Austin.  The  flow  diagram 
of  the  program  used  in  the  computation  of  the  ROC  curves  is  listed 
in  Appendix  B. 

B.  Target  Classification  for  Various  S/N  Ratios 

In  this  section,  the  performance  of  the  likelihood  ratio 
processor  in  classifying  targets  using  the  s/N  ratio  of  the  messages 
from  the  target  classes  as  a variable  parameter  was  measured.  The 
transmission  mode  used  was  the  sinusoidal  pulse  described  in  Chapter  III. 
The  rms  of  the  amplitude  of  the  noise  to  be  added  to  the  echoes  from  a 
target  in  order  to  lower  the  S/N  ratio  of  the  messages  from  the  target 
class  to  a desired  value  was  specified  by  Eq.  (4.2).  The  knowledge  of 
epoch  was  retained  during  the  addition  of  the  noise  to  the  echoes  so 


-;at  the  techniques  discussed  in  Chapter  IV,  Section  A,  and  the  test 

n 

.A 


V,  s®  + s®  - V,  - v„  s;: 

Im  In  2m  2n  Im  In  2m  2n 


could  be  used  in 


■ flcation  of  the  targets. 
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The  following  three  pairs  of  targets  were  used:  the  7 in. 

diam  and  5 in-  diam  solid  aluminum  spheres,  the  5 in.  diam  and  3 in. 
diam  solid  aluminum  spheres,  £ind  the  5 in.  diam  styrofoam  sphere  and 
the  5 in.  diam  hollow  aluminum  sphere  (described  in  Chapter  III). 

Four  ROC  curves  for  each  target  pair  were  computed,  each  curve  corre- 
sponding to  a different  S/N  ratio  of  the  messages  from  each  of  the  two 
target  classes.  The  four  values  for  the  S/n  ratio  were  6 dB,  3-5  dB, 

0 dB,  and  -6  dB.  The  ROC  curves  are  plotted  in  Figs.  4.11  through  4.13. 


Txie  area  under  each  ROC  curve  along  with  the  values 


'"A  - "b' 


, jr , and  — are  listed  in  Tables  4.6  through  4.8. 

p o 


— multiplied  by  100  gives  the  percent  overlap  of  the  hyperspheres 
^T 

generated  by  the  vectors  and  in  a 2TW' -dimensional  vector  space. 
For  each  target  pair,  the  percent  overlaps  for  the  four  S/N  ratios  are 
listed  in  Table  4.9- 


C.  Target  Classification  Using  One  Quadrature  Component 

In  this  section  the  ability  of  a likelihood  ratio  processor 
to  classify  messages  of  known  epoch  using  only  one  quadrature  component 
is  measured.  The  two  targets  used  were  the  in.  diam  and  5 in.  diam 
solid  aluminum  spheres.  The  waveform  used  in  transmission  was  the 
sinusoidal  pulse  described  in  Chapter  III.  From  Eq.  (2.69) j the  equation 
for  i(v. ) using  only  the  in-phase  components  and  is 
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Q 

- ^ a ^ fv  - V-i  S®  - 1/2  + 1/2  S®  S®  1 
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PROBABILITY  OF  CORRECTLY  CLASSIFYING  TARGET 


0 0.2  0.4  0.6  0.8 

PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  S/N  RATIO 

S/N  * 6 dB 

S/N  = 3.5  dB 

S/N  = 0 dB 

S/N  = -6  dB 


FIGURE  4.11 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  VARIOUS  S/N  RATIOS 
TARGET  A * 7 in.  diam  SOLID  ALUMINUM  SPHERE 
TARGET  B = 5 in.  diom  SOLID  ALUMINUM  SPHERE 


TABLE  4.6 

TARGET  A = 7 IN.  DIAM  SOLID  ALUMINUM  SPHERE 
TARGET  B = 5 IN.  DIAM  SOLID  ALUMINUM  SPHERE 


r -r 
A B 


VS 


Area  under 
ROC  Curve 


S/'N  = 6 dB 

S'N  = 3.5  dB 

S ;'N  = 0 dB 

S/N  = -6  dB 

0.03 

0.02 

0.009 

0.004 

0.09 

0.08 

0.06 

0.04 

1.65 

1.46 

1.29 

1.13 

0.899 

0.844 

0.760 

0.644 
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CODE  FOR  S/N  RATIO: 

S/N  = 6 dB 

S/N  = 3.5  dB 

S/N  = OdB 

S/N  = -6  dB 


FIGURE  4.12 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
^ ' FOR  VARIOUS  S/N  RATIOS 

{ TARGET  A s 5 in.  diam  SOLID  ALUMINUM  SPHERE 

TARGET  B = 3 in.  diam  SOLID  ALUMINUM  SPHERE 


I 

j 


TARGET  A = 5 IN.  DIAM  SOLID  ALUMINUM  SPHERE 
TARGET  B = 3 IN.  DIAM  SOLID  ALUMINUM  SPHERE 


S/N  = 6 dB  S/N  = 3.5  dB  S/N  = 0 dB  S 'N  = -6  dB 


A B 


Area  under 
ROC  Curve 
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PROBABILITY  OF  CORRECTLY  CLASSIFYING  TARGET  A 


PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  S/N  RATIO: 

S/N  = 6 dB 

S/N  = 3.5  dB 

S/N  = OdB 

S/N  = -6dB 


FIGURE  4.13 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  VARIOUS  S/N  RATIOS 
TARGET  A « 5 in.  diom  STYROFOAM  SPHERE 
TARGET  B = 5 in.  diam  HOLUDW  ALUMINUM  SPHERE 
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Tho  equation  for  using  only  the  out-of-phase  components 


and  ^2  is 


,2  ^ w; 


£(^2)  = e ^ 


. (^.7) 


From  Eq.  (2.69),  Eq.  (4.6),  ajid  Eq.  (4.7),  it  follows  that 


i(^  = i(^^)  , or 

log  (^(V^]  = log  + log  {>(%)) 


(^.8) 

(i^.9) 


If  log  (i(^^))  and  log('i(f2))  are  of  opposite  sign,  the  sign  of 
log^i  C\f)j  will  he  determined  by  the  larger  of  the  two  numbers,  | log ^^(V^) j ) , 
or  ilog^i(?2))|  • 

In  lh.g.  4.l4  the  in-phase  and  the  out-of-phase  components  of 
an  echo  from  a 7 in.  diam  solid  aluminum  sphere  and  the  in-phase  and 
the  out-of-phase  components  of  an  echo  from  a 5 in.  diam  solid  aluminum 
sphere  are  plotted.  In  Fig.  4.14  the  S/n  ratio  of  the  echo  from  the 
7 in.  diam  solid  aluminum  sphere  was  59-8  dB,  and  the  S/N  ratio  of  the 


echo  from  the  5 in.  diam  solid  aluminum  sphere  was  38.8  dB.  In  Fig.  4.15 
the  in-phase  and  the  out-of-phase  components  of  a message  from  each  of 
the  two  target  classes  are  plotted.  The  S/N  ratio  of  the  messages  was 
6 dB.  The  ROC  curves  in  Fig.  4.l6  and  in  Fig.  4.17  were  plotted  using 


r 
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0.37  msec 
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! 

OUT-OF-PHASE  COMPONENT  OF  ECHO 
FROM  5 in.  diam  SOLID  ALUMINUM  SPHERE 
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IN-PHASE  COMPONENT  OF  ECHO 
FROM  5 in.  diam  SOLID  ALUMINUM  SPHERE 
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OUT-OF-PHASE  COMPONENT  OF  ECHO 
FROM  7 in.  diam  SOLID  ALUMINUM  SPHERE 
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IN-PHASE  COMPONENT  OF  ECHO 
FROM  7 in.  diam  SOLID  ALUMINUM  SPHERE 


FIGURE  4.14 

QUADRATURE  COMPONENTS  FOR  SINUSOIDAL  PULSE 

PULSE  LENGTH  = 0.37  msec  EPOCH  = = 8/300  msec 

= 0 kHz  W = 20  kHz  PEAK  |AMPLITUDE|  = 0.5  in. 
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OUT-OF-PHASE  COMPONENT  OF  ECHO 
FROM  5 in.  diam  SOLID  ALUMINUM  SPHERE  PLUS  NOISE 
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IN-PHASE  COMPONENT  OF  ECHO 
FROM  5 in.  dio-m  SOLID  ALUMINUM  SPHERE  PLUS  NOISE 
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OUT-OF-PHASE  COMPONENT  OF  ECHO 
FROM  7 in.  diam  SOLID  ALUMINUM  SPHERE  PLUS  NOISE 
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IN-PHASE  COMPONENT  OF  ECHO 
FROM  7 in.  diam  SOLID  ALUMINUM  SPHERE  PLUS  NOISE 


FIGURE  4.15 

QUADRATURE  COMPONENTS  OF  MESSAGES  FOR  A SINUSOIDAL  PULSE 
POSITION  OF  0.37  msec  TIME  INTERVAL  INDICATES  LOCATION  OF  EMBEDDED  COMPONENT 
PULSE  LENGTH  = 0.37  msec  At  = 8/300  msec  = 0 kHz 

W = 20  kHz  S/N  = 6 dB  PEAK  |AMPLITUDE|  = 0.5  in. 
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FIGURE  4.16 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
COMPUTED  FROM  ONE  QUADRATURE  COMPONENT 
AND  BOTH  QUADRATURE  COMPONENTS 
TARGET  A = 7 in.  diam  SOLID  ALUMINUM  SPHERE 
TARGET  B = 5 in.  diam  SOLID  ALUMINUM  SPHERE 
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S/N  = 0 dB 


PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  QUADRATURE  COMPONENTS; 

IN-PHASc  PLUS  OUT-OF-PHASE  COMPONENTS 

IN-PHASE  COMPONENT 

OUT-OF-PHASE  COMPONENT 


FIGURE  4.17 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
COMPUTED  FROM  ONE  QUADRATURE  COMPONENT 
AND  BOTH  QUADRATURE  COMPONENTS 
TARGET  A 7 in.  diam  SOLID  ALUMINUM  SPHERE 
TARGET  B - 5 in.  diam  SOLID  ALUMINUM  SPHERE 
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the  test  functions 
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respectively  and  the  techniques  discussed  in  Chapter  IV,  Section  A. 

The  test  functions  are  monotonically  increasing  functions  of 

and  respectively.  The  average  s/N  ratio  of  the  messages 

from  each  target  class  in  Fig.  4.l6  was  6 dB,  and  the  average  s/N 
ratio  of  the  messages  from  each  target  class  in  Fig.  4.17  was  0 dB. 

The  area  under  each  ROC  curve  in  Fig.  4.l6  and  Fig.  4.17  is  listed 
in  Table  4.10. 

From  Fig.  4.l6  and  Fig.  4.17,  the  test  function  that  used 
the  in-phase  components  gave  better  performance  than  the  test  function 
that  used  the  out-of-phase  components.  No  general  statement  can  be 
made  from  this  result  since  it  is  highly  dependent  on  the  transmission 
mode,  epoch,  and  the  target  pair. 

The  best  performance  in  classification  was  from  the  test 
function  that  used  both  quadrature  components.  This  result  was  pre- 
dicted in  Eq.  (4.9). 
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D.  Target  Classification  Usinfi  a Generalized  Likelihood 
Ratio  Processor 

In  this  section  the  performajrice  of  the  generalized  likelihood 
ratio  processor  ^Eq.  (2.6T))  is  compared  with  the  performance  of  the 
likelihood  ratio  processor  ^Eq.  (2.69)j  in  the  classification  of  tar- 
gets. The  two  targets  used  were  the  7 in.  diam  and  5 in.  diam  solid 
aluminum  spheres.  The  average  s/n  ratio  of  the  messages  from  each 
target  class  was  6 dB.  The  starting  phase  of  the  messages  was  varied 
in  intervals  of  it/2  radians  over  a full  cycle  of  the  carrier.  The 
number  of  values  of  epoch  for  each  message  was  thereby  increased  from 

( B^W^  I ^ 

1 to  4.  The  test  function  +v/  -\\  +v/  , which  is  a 

monotonically  increasing  function  of  described  by  Eq.  (2.67)  where 


TW- 


w; 


V,  s®  + V s® 
Ip  Iq  2p  2q 


etc.,  was  calculated  for  each  value 


of  epoch.  Next,  an  average  value  for  the  test  functions  was  computed 
for  each  message.  In  Fig.  4.l8  an  ROC  curve  was  plotted  using  the  aver- 
age values  for  the  test  functions.  A second  ROC  curve  was  computed  for 
the  same  target  pair  and  was  plotted  on  the  same  graph.  The  second  ROC 
curve  was  computed  from  messages  of  known  epoch  using  the  test  function 


p,q=l 


which  is  a monotonically 


increasing  fimcticn  of  described  by  Eq.  (2.69). 

The  degradation  in  performance  due  to  the  loss  of  epoch  can 
be  estimated  frpjj^^comparlson  of  the  two  curves.  For  the  likelihood 
ratio  processor,  the  probability  of  correctly  classifying  target  A 
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PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


i KNCWN  EPOCH 

/ AVERAGE  OP  VARIED  EPOCH 

I 


FIGURC  4.18 

^ ^ RECEIVER  OfERATING  CHARACTCttSTIC  CURVES 

I FOR  KNOWN  EPOCH  AND  AVERAGE  OF  VARIED  EPOCH 

^ . TARGET  A • 7 in.  diom  SOLID  ALUMINUM  SPHERE 

i TARGET  B < 5 in.  diom  SOLID  ALUMMUM  SPHERE 


was  86^  when  the  prohahility  of  incorrectly  classifying  target  B was 
20^,  and  for  the  generalized  likelihood  ratio  processor,  the  probability 
of  correctly  classifying  target  A was  6o^  when  the  probability  of 
incorrectly  classifying  target  B was  20^.  From  Fig.  4.11  the  degra- 
dation in  the  performance  of  classification  because  of  loss  of  epoch 
of  the  messages  is  approximately  equal  to  the  degradation  in  the 
performance  due  to  a 6 dB  decrease  in  the  S/N  ratio  of  the  messages. 


CHAPTER  V 


TARGET  CLASSIFICATION  USING  L'TJEAR  FW  TRANSIISSION 


A.  Classification  of  Similar  Targets 
1 . Epoch 

The  investigation  concerning  the  classification  of  targets 
using  a linear  FM  transmission  is  presented  in  this  chapter.  A plot 
of  the  linear  M transmission  is  shown  in  Fig.  3.2.  The  echoes  from 
the  linear  EM  transmission  were  passed  through  a 5 msec  gate.  A 
plot  of  a gated  echo  from  each  of  the  5 targets,  the  7 in.  diam,  the 
5 in.  diam,  and  the  3 in.  diam  solid  aluminum  spheres,  the  5 in.  diam 
styrofoam  sphere,  and  the  5 in.  diam  hollow  aluminum  sphere,  is  shown 
in  Fig.  5.k. 


The  epoch  for  a gated  echo  was  detennined  by  the  time,  t^, 

F F 

at  which  Z ( t)  had  a maximum  value,  where  Z (t)  is  the  correlation 

between  a gated  echo  ^ and  a transmitted  signal  The  epoch  of  an 

echo  was  defined  to  be  e , where  e = t . The  epoch  for  an  echo  from 

o o o 

each  target  is  indicated  in  Fig.  J>.h, 

The  correlation  between  a gated  echo  ^ and  a transmitted 
^T 

waveform  S for  N,.  + 1 number  of  lags  is  described  as  follows: 

Jj 


and 


N. 


a+T 


0=1 


(5.1) 


T = 0,  1,  2, 
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ft.  . 


target.  In  Table  5.1  the  following  three  qmntities  for  each  target 
are  tabulated: 


95 


(2)  D = 


1 1/2 


T ^ 


L “ J=1  \<jV^  > 


172 


- M 


, and 


(5.2) 


(5.5) 


(5)  ^ X 100 


where 


= l400  = peak  amplitude  of  each  echo  from  the  target, 

< = rms  of  the  amplitude  of  echo  from  the 

target,  and 

N = number  of  echoes  from  the  target . 

i!i 

The  quantity  D/^  is  small  for  each  target. 

2.  Addition  of  Gaussian  Noise  to  Echoes 

Gaussian  noise  frou  the  random  noise  generator  was  added 
to  the  echoes  from  each  target  in  order  to  meike  the  simulation  of 
target  classification  more  realistic  and  to  increase  the  number  of 
messages  (echo-noise  combinations)  frcan  each  target  class.  The 
bandwidth  of  the  Gaussian  noise  was  20  kHz  with  a center  frequency 
of  65  kHz . The  rms  of  the  amplitude  of  the  noise  from  the  random 
noise  generator  was  specified  as  follows: 


< N > 


2 . 1/2  _ P, 


R, 


(5.U) 


SN 
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where 


2 1/2 

< N > ' = rms  of  the  amplitude  of  the  noise  from  the  random 

noise  generator, 

= 1^400  = peak  amplitude  of  each  echo  from  the  target,  and 
Rgjj  = desired  S/H  ratio  of  the  echoes  from  the  target , 

The  rms  of  the  amplitude  of  the  ambient  noise  and  of  the  reverberation 
present  in  each  echo  was  considered  to  be  negligible  compared  to  the 
peak  amplitude  of  each  echo.  In  Fig.  5.1  a message  with  a S/N  ratio 
of  -6  dB  from  each  target  class  is  plotted.  In  Fig.  5.1  the  messages 
of  pulse  length  5.2  msec  are  located  between  the  two  vertical  lines. 

A comparison  of  the  plots  in  Fig.  5.^  and  Fig.  5.1  shows  that  it  is 
very  difficult  to  classify  visually  the  echoes  that  are  embedded  in 
Gaussian  noise. 

For  the  linear  M transmission  the  number  of  echoes 
available  from  each  target  varied  from  16  to  I8.  After  the  addition 
of  noise  to  the  echoes,  the  number  of  messages  from  each  target  class 
varied  from  128  to  l44.  No  noise  signal  was  duplicated  throughout 
the  procedure  involving  the  addition  of  noise  to  the  echoes . 

In  obtaining  the  quadrature  components  of  the  messages, 
the  center  frequency  of  the  messages  was  redefined  as  f^  = 75  kHz, 
and  the  bandwidth  was  redefined  as  W = 40  kHz.  A sampling  parameter 
of  K = 1 was  used,  which  resulted  in  a sampling  rate  W'  of  75  WJz. 

The  likelihood  ratio  calculation  described  by  Eq.  (2.69) 


was  simplified  by  approximating  the  matrix 


as  a diagonal  matrix; 


Xfi*. 
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namely. 


(■5)  . 


(5.5) 


where 


WA  w- 

a.?  = a ^ 
ij 


The  approximation 


W ■ 
(•5)-( 


2 

< N > 


.) 


^ijj  necessary  since  the  digital 


processing  involving  a matrix  of  the  size  TW'  by  TW' , where  TW'  is 
2k0,  was  not  feasible.  Justification  for  the  use  of  a diagonal 
matrix  was  that  an  ROC  curve  computed  using  a diagonal  matrix  and 
messages  for  a sinusoidal  pulse  did  not  differ  measurably  from  an  ROC 
curve  computed  using  a nondiagonal  matrix  and  the  same  messages.  The 

1 


approximation 


o - S . . I was  not  needed  in  Chapter  TV, 

V ^ J / 


<<  N > 

since  the  size  of  the  matrix  required  for  the  digital  processing  of 
the  messages  for  the  sinusoidal  pulse  did  not  present  any  computa- 
tional difficulties. 

Upon  the  substitution  of  Eq.  (5.5)  into  Eq.  (2.69),  the 
likelihood  ratio  used  in  the  classification  of  the  messages  for  a 
linear  M transmission  is  described  as  follows: 


Tw: 


'2W^  S (<  ^2  >)  ""  ^2p4p  - \p^lp  ■ ^2p^2p 


p=l 


- - Sp  Sp  + s®  S®  + Sp  s® 

Ip  Ip  2p  2p  Ip  Ip  2p  2p  , 


(5.6) 


«. 


where 
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where  TW^  is  the  number  of  samples  of  each  quadrature  component . 

3 . ROC  Curves 

Four  ROC  curves  are  plotted  in  Fig.  5.2  for  the  linear  M 
transmission.  These  curves  were  computed  using  the  test  function 

% r I 

E " ''pp=®p  - - ''pp4pj  . 

p=i 

This  test  function  is  a monotonically  increasing  function  of 
where  is  described  by  Eq.  (5.6).  In  Fig.  5.2,  target  A for  each 

of  the  4 ROC  curves  was  the  7 in.  diam  solid  aluminum  sphere,  and 
the  average  S/M  ratio  of  the  messages  frcan  each  target  class  was 
-6  dB. 

A comparison  of  the  ROC  curves  in  Fig.  5.2  with  those  in 
Fig.  4.6,  which  were  computed  for  the  sinusoidal  pulse,  reveals 
that  the  ROC  curves  for  the  different  spheres,  B,  are  in  the  same 
order.  But  each  curve  in  Fig.  5.2  indicates  a better  performance 
in  classification  than  the  corresponding  curve  in  Fig.  4.6,  even 
though  the  average  S/W  ratio  of  the  messages  used  in  the  computation 
of  the  ROC  curves  in  Fig.  4.6  was  0 dB. 
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PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  TARGET  B 

5 in.  diom  SOLID  ALUMINUM  SPHERE 

3 in.  diom  SOLID  ALUMINUM  SPHERE 

5 in.  dkjm  STYROFOAM  SPHERE 

5 in.  diom  HOLLOW  ALUMINUM  SPHERE 


FIGURE  5.2 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  LINEAR  FM  TRANSMISSION 
TARGET  A * 7 in.  diom  SOLID  ALUMINUM  SPHERE 
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B.  Conditions  Necessary  for  EqulVEilent  Performance  in  Classification 
Using  Either  of  the  Two  Transmissions 

The  reason  for  the  better  results  Indicated  by  Fig.  5.2  is 
that  an  echo  from  the  linear  M transmission  contained  more  energy 
than  an  echo  from  the  sinusoidal  pulse.  Since  the  energy  in  an  echo 
is  proportional  to  the  pulse  length  of  the  echo  times  the  peak  ampli- 
tude squared  and  since  every  echo  had  the  same  peak  amplitude,  the 
greater  energy  of  an  echo  from  the  linear  M transmission  was  the 
result  of  the  longer  pulse  length  of  the  linear  M transmission. 

The  reciprocal  of  the  pulse  length  of  the  sinusoidal  pulse 
was  equal  to  the  bandwidth  of  the  linear  M transmission.  Therefore, 
from  simple  energy  considerations  and  by  ignoring  waveform,  one  would 
expect  equivalent  perfomances  in  classification  if  the  following 
equality  were  satisfied: 


where 


E = energy  in  an  echo, 

2 

< N > = mean  square  of  the  amplitude  of  the  added  Gaussian 
noise, 

sp  denotes  the  sinusoidal  pulse,  and 
PW  denotes  the  linear  EM  transmission. 

Since  the  energy  in  an  echo  is  proportional  to  the  pulse 
length  of  the  echo  times  the  peak  amplitude  squared,  Eq.  (5.7)  can 


be  written  as 


\<  \<  >L 


'■  ilL 


1 


where 


T = pulse  length  of  an  echo  and 
= l400  = peEik.  amplitude  of  an  echo. 


Upon  taking  the  square  root  of  Eq.  (5.8),  the  result  is 

/ \ / \ 


^ 1/2 


^<N^>  7 V 

' 'sp  ' 


tV^Pa  \ 

N > /v 


(5.9) 


2 

< N > 


the  echo-noise  combinations  for  the  sinusoidal  pulse,  the  result  is 


After  solving  Eq.  (5.9)  for  j — jTpj  > ratio  of 


2 

K N > 


\ _ /^m\  / 

2 I I T I I „ 1 

/s  V ^ > 


, or,  in  dB,  (5.10) 


20  log 


1a)  “ ' “ . 20  log 

> /sp  V°7  > L 


(5.11) 


From  Eq.  (5.H),  20  log 


V 2 V 

\<  N > 


dB  is  equal  to  5.^  dB 


when  20  log 


V<  > 


-r-fir  I dB  = -6  dB,  T = 0.57  msec,  and  T_,  = 3.2  msec. 
1/21  ' sp  M 


In  Fig.  5.3,  ^ ROC  curves  were  plotted  for  the  sinusoidal  pulse.  The 
s/s  ratio  of  the  messages  used  in  the  computation  of  the  ROC  curves 
was  3.4  dB.  The  performances  indicated  by  the  ROC  curves  in  Fig.  5.3 
duplicated  the  performances  indicated  by  the  ROC  curves  in  Fig.  5.2. 
The  area  under  each  ROC  curve  in  Figs.  5.2  and  5.3  is  listed  in 


Table  5.2. 
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PROBABILITY  OF  INCORRECTLY  CLASSIFYING  TARGET  B 


CODE  FOR  TARGET  B 

5 in.  diam  SOLID  ALUMINUM  SPHERE 
3 in.  diam  SOLID  ALUMINUM  SPHERE 
5 in.  diam  STYROFOAM  SPHERE 
5 in.  diam  HOLLOW  ALUMINUM  SPHERE 


FIGURE  5.3 

RECEIVER  OPERATING  CHARACTERISTIC  CURVES 
FOR  SINUSOIDAL  PULSE  TRANSMISSION 
TARGET  A * 7 in.  diam  SOLID  ALUMINUM  SPHERE 
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TABLE  5.2 

AREA  UNDER  ROC  CURVES  IN  FIGURES  5.2  AND  5.3 
TARGET  A = 7 IN.  DIAM  SOLID  ALUMINUM  SPHERE 


Target  B 

5 in.  diam 
Solid  Aluminuir 
Sphere 

3 in.  diatn 
Solid  Aluminum 
Sphere 

5 in . diam 
Styrofoam 
Sphere 

5 in . diam 
Hollow  Aluminum 
Sphere 

Figure  5*2 

0.838 

0.861 

0.876 

0.864 

Figure  5-3 

0.842 

0.862 

0.900 

0.865 

I 
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CHAPTER  VI 


CONCLUSIONS 

The  objective  of  this  study  was  the  quantitative 
determination  of  the  applicability  of  a likelihood  ratio  processor  in 
the  classification  of  sonar  returns  from  targets . The  ability  of  a 
likelihood  ratio  processor  to  classify  simple,  similar  targets  was 
measured  under  highly  controlled  conditions . The  5 targets  used  were 
the  7 in.  diam,  5 in.  diam,  and  3 in.  diam  solid  aluminum  spheres,  the 
5 in.  diam  styrofoam  sphere,  and  the  5 in.  diam  hollow  aluminum 
sphere,  as  described  in  Chapter  III.  The  two  transmissions  used  were: 

(1)  a sinusoidal  pulse  with  a frequency  of  70  kHz  and  a 
pulse  length  of  80  psec,  and 

(2)  a linear  M transmission  of  bandwidth  12.5  kHz  with  a 
center  frequency  of  65.75  kHz  and  a pulse  length  of  3 msec. 

The  echoes  from  the  targets  for  both  modes  of  transmission 
were  scaled  so  that  each  echo  had  the  same  peak  amplitude.  The 
purpose  of  this  scaling  was  to  be  able  to  base  any  future  discrimi- 
nation solely  on  the  waveform  of  the  echoes;  namely,  echo  strength 
was  not  to  be  a factor  in  classification.  For  the  echoes  from  the 
sinusoidal  p’Jilse,  the  epoch  was  determined  by  a threshold  technique; 
for  the  echoes  from  the  linear  FM  transmission,  the  epoch  was 
determined  by  the  time  at  which  a correlation  function  obtained  its 
maximum  value. 
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In  order  to  simulate  a more  realistic  problem  in  target 
classification  and  to  increase  the  number  of  messages  (echo-noise 
combinations)  from  each  target  class,  Gaussian  noise  from  a random 
noise  generator  was  added  to  the  echoes.  Upon  the  addition  of 
noise  to  the  echoes,  the  knowledge  of  epoch  was  retained.  It  was 
possible  to  select  the  S/N  ratio  of  the  messages  by  specifying  the 
rms  of  the  amplitude  of  the  noise  from  the  random  noise  generator. 

In  classifying  messages  from  two  target  classes,  the  performance  of 
the  likelihood  ratio  processor  was  measured  by  means  of  ROC  curves , 
The  number  of  messages  used  in  computing  an  ROC  curve  was  about  15O 
for  each  target  class . 

For  the  sinusoidal  pulse,  the  likelihood  ratio  processor 
performed  satisfactorily  in  classifying  the  messages  for  a majority 
of  the  target  pairs  when  the  S/N  ratio  of  the  messages  was  0 dB  or 
larger.  A performance  was  considered  satisfactory  if  the  area  under 
an  ROC  curve  was  0.75  or  greater,  since  the  area  under  an  ROC  cvirve 
times  100  is  equal  to  the  expected  percentage  of  correct  decisions. 
For  a S/N  ratio  of  0 dB,  it  was  virtually  impossible  to  classify  or 
to  detect  visually  the  echoes  embedded  in  noise . 

An  alternate  method  of  measuring  the  ability  to  classify 
messages  utilized  a geometric  representation  of  the  messages  frcan 
each  target  class  in  a multidimensional  vector  space.  For  a given 
target  A and  target  B,  the  degree  of  overlap  of  the  vector  space 
containing  the  messages  from  target  class  A with  the  vector  space 
containing  the  messages  from  target  class  B was  indicative  of  the 


i 


r 


1 
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performance  of  the  likelihood  ratio  processor.  Namely,  the  smaller 
the  overlap,  the  larger  was  the  area  under  the  ROC  curve. 

For  the  sinusoidal  pulse,  the  effects  of  variations  in  the 
ratio  of  the  messages  on  the  performance  of  the  likelihood  ratio 
processor  were  measured.  The  S/iJ  ratio  of  the  messages  was  varied 
from  6 dB,  3.5  dB,  O dB,  to  -6  dB.  The  S/ta  ratio  was  a sensitive 
parameter  that  affected  the  performance  of  the  likelihood  ratio  pro- 
cessor. The  larger  the  S/Is  ratio,  the  larger  was  the  area  under  the 
ROC  curve.  As  was  expected,  the  degree  of  overlap  of  the  vector 
spaces  generated  by  the  messages  was  also  dependent  on  the  S/^  ratio. 

The  larger  the  ratio,  the  smaller  was  the  overlap  of  the  vector 
spaces . 

For  the  sinusoidal  pulse,  the  effect  of  a loss  in  epoch  on 
the  ability  to  classify  the  messages  was  measured.  The  loss  it.  epoch 
occurred  when  the  starting  phase  of  the  messages  was  varied  in 
intervals  of  «/2  radians  over  a full  cycle  of  the  carrier.  A test 
function  that  was  a monotonically  increasing  function  of  a generalized 
likelihood  ratio  was  calculated  for  each  value  of  epoch,  and  for  each 
message  an  average  test  function  was  computed.  The  loss  in  epoch  had. 
the  same  effect  on  the  ability  to  classify  the  messages  as  a 6 dB 
drop  in  the  S/lil  ratio  of  the  messages . The  degradation  in  performance 
was  due  to  a loss  in  coherency  between  the  messages  and  their  re- 
spective reference  signals. 

For  the  sinusoidal  pulse,  it  was  possible  to  classify 
messages  with  a likelihood  ratio  processor  that  vised  only  in-phase 
components  or  out-of- phase  components.  For  target  A equal  to  the 


7 in.  diam  solid  aluminum  sphere  and  target  B equal  to  the  5 in.  diam 
solid  aluminum  sphere,  the  in-phase  components  resulted  in  better 
performance  in  classification  than  the  out-of- phase  components.  But 
this  result  is  limited,  in  that  it  is  dependent  on  the  arbitrary 
manner  in  which  the  quadrature  components  were  defined.  Better  per- 
formance in  classification  was  obtained  with  both  quadrature  com- 
ponents than  with  only  one  quadrature  component,  since  the  knowledge 
of  both  components  gave  the  likelihood  ratio  processor  additional 
Information  on  which  to  base  its  classification. 

The  performance  obtained  with  the  likelihood  ratio 
processor  as  indicated  by  ROC  curves  for  the  linear  M transmission 
when  the  S/iT  ratio  of  the  messages  was  -6  dB  was  better  than  that 
obtained  for  the  sinusoidal  pulse  when  the  ratio  of  the  messages 
was  0 dB.  The  reason  for  the  better  performance  was  that  an  echo 
from  a linear  M transmission  contained  more  energy  than  an  echo 
from  a sinusoidal  pulse. 

Since  the  reciprocal  of  the  pulse  length  of  the 
sinusoidal  pulse  was  equal  to  the  bandwidth  of  the  linear  PM,  the 
two  transmissions  had  equivalent  capabilities  in  resolving  the 
separate  components  of  an  echo.  Therefore  from  simple  energy 
considerations  and  ignoring  waveform  it  was  plausible  that  equivalent 
performances  could  be  obtained  with  either  of  the  two  transmissions 
if  the  following  equation  were  true: 


108 


> 

I 


it 

r 4 


\<  4p  ' V >)m 


(6.1) 


where 


T = pulse  length  of  an  echo, 

= 1^+00  = peak  amplitude  of  an  echo, 

2 

E is  proportional  to  TP^  where  E is  the  energy  in  an  echo, 

2 

< N > = mean  square  of  amplitude  of  noise  from  the  random 
noise  generator, 

sp  denotes  a sinusoidal  pulse,  and 
M denotes  a linear  M transmission. 

The  performances  achieved  with  the  linear  transmission 
when  the  ratio  of  the  messages  was  -6  dB  were  also  achieved  with 
the  sinusoidal  pulse  when  the  S/n  ratio  of  the  messages  was  5.4  dB 


The  value  5.4  dB  was  obtained  by  solving  Eq.  (6.1)  for 


the  ratio  of  the  messages  for  the  sinusoidal  pulse,  when 
was  equal  to  -6  dB. 


V"  ^ /sp 


f "a  \ 
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From  Eq.  (6.l)  it  was  asserted  that  the  performance  of  the 

E 


likelihood  ratio  processor  is  dependent  on  the  value  of 


namely,  the  larger 


E 


2 ^ 
< N > 


2 

< N > 


is,  the  better  is  the  performance.  The 


dependence  of  the  performance  on 


E 


implies  that  a linear  M 


< > 

has  more  potential  than  a sinusoidal  pulse  for  improvement  in  the 
perfcmance  of  a likelihood  ratio  processor,  since  the  energy  in  a 
sinusoidal  pulse  is  restricted  by  the  effects  of  cavitation. 
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RANDOM  NOISE  GENERATOR 


The  random  noise  generator  used  in  obtaining  samples  of 
Gaussian  noise  was  the  type  1390-B  rajidom  noise  generator  manufactured 
by  General  Radio  Company.  The  type  1J90-B  random  noise  generator 
uses  a gas  discharge  tube  as  its  noise  source.  Samples  of  the  ampli- 
tude of  the  output  waveform  are  characterized  by  a Gaussian  distribution. 

28 

Two  of  the  pertinent  specifications  on  the  type  I39O-B 
random  noise  generator  are: 

(1)  frequency  range:  5 Hz  to  5 MHz  and 

(2)  spectrum-level  uniformity:  within  ±3  dB  from  20  Hz 

to  500  kHz. 


28 


Operating  Instructions,  Type  1390-B  Random  Noise  Generator, 
General  Radio  Company,  West  Concord,  Mass.,  February  196^. 
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APPENDIX  B 

FLOW  DIAGRAM  OF  COMPUTER  PROGRAM 
USED  IN  CALCULATION  OF  ROC  CURVE 
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APPENDIX  B CONTINUED 


1K 


APPENDIX  B CONTINUED 


Ak  = INCREMENT  OF 
THRESHOLD 


BIBLIOGRAPHY 


Bendat,  J.  S.,  ajrid  A.  G,  Pier  sol.  Measurement  and  Analysis  of 
Random  Data.  New  York:  John  Wiley  and  Sons,  Inc.,  19^6. 

Blackman,  R.  B.,  and  J.  W.  Tukey.  The  Measurement  of  Power  Spectra. 

New  York:  Dover  Publications,  Inc.,  19 9^- 

Dierks,  K.  Jerome,  and  R.  Hickiing.  "Echoes  from  Hollow  Aluminusi 
Spheres  in  Water."  Journal  of  the  Acoustical  Society  of 
America  No.  2 (196?),  pp.  580-393. 

Grace,  0.  D.,  and  S.  P.  Pitt.  "Quadrature  Sampling  of  High  Frequency 
Waveforms."  Letter  to  the  Editor  of  the  Journal  of  the  Acoustical 
Society  of  America  No.  5 (1968),  pp. 

Green,  David  M.,  and  John  A.  Swets.  Signal  Detection  Theory  and 
Psychophysics.  New  York:  John  Wiley  and  Sons,  Inc.,  19^(3. 

Hancock,  J.  C.,  and  P.  A.  Wintz.  Signal  Detection  Theory. 

New  York:  McGraw-Hill  Book  Company,  1966. 

Hickiing,  Robert.  "Analysis  of  Echoes  from  a Hollow  Metallic  Sphere 
in  Water."  Journal  of  the  Acoustical  Society  of  America  36, 

No.  6 (196^+),  pp.  112i+-1137. 

Hickiing,  Robert.  "Analysis  of  Echoes  from  a Solid  Elastic  Sphere 
in  Water."  Journal  of  the  Acoustical  Society  of  America  3^> 

No.  10  U962),  pp.  1582-1592. 

Hickiing,  Robert,  and  N.  M.  Wang.  "Scattering  of  Sound  by  a Rigid 
Movable  Sphere."  Journal  of  the  Acoustical  Society  of  America 
No.  2 (1966),  pp.  276-279. 

Middleton,  David.  An  Introduction  to  Statistical  Communication 
Theory.  New  York:  McGraw-Hill  Book  Company,  i960 . 

Middleton,  David.  Topics  in  Communication  Theory.  New  York: 
McGraw-Hill  Book  Company,  1965. 

Papoulis,  Athanasious.  "Limits  on  Band-limited  Signals."  Proceedings 
of  the  IEEE  55,  No.  10  (1967),  pp.  1677-1686. 

Schooneveld,  C.  Van.  "Some  Remarks  on  Sampling  Methods  for  a Band 

Pass  Signal."  Signal  Processing,  NATO  Advanced  Study  Institute, 
Grenoble  (196^),  pp.  ^10-421. 

Sebestyen,  G.  S.  Decision-making  Processes  in  Pattern  Recognition. 

New  York:  The  Macmillan  Company,  I962. 

115 


June  1969 


DISTRIBUTION  LIST  FOR 
ARL-TM-69-11 

CONTRACT  N0002lf-69-C-1129 
Project  Serial  No.  SF  11121100,  Task  O8515 


Commander,  Naval  Ship  Systems  Command 
Department  of  the  Navy 
Washington,-  D.  C.  20j60 
Attn:  SHIPS  OOVIC 

Commander,  Naval  Ship  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20j60 
Attn:  SHIPS  00V2 

Commander,  Naval  Ship  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20j60 
Att : PMS  385 

Commander,  Naval  Ship  Systems  Command 
Department  of  the  Navy 
Washington , D . C . 203^0 

Attn : PMS  386 

Commander,  Naval  Ship  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  2036O 
Attn : PMS  38? 

Commander,  Naval  Ship  Systems  Command 
Department  of  the  Na-vy 
Washington,  D.  C.  2036O 
Attn:  Library 

Chief  of  Naval  Research 
Department  of  the  Navy 
Washington,  D.  C.  2036O 
Attn:  Code  466 


116 


117 


June  1969 


DISTRIBUTION  LIST  FOR 
ARL-TM-69-11 

CONTRACT  NOOO24-69-C-II29  (CONT'D) 


Commander 

Naval  Undersea  Research  and  Development  Center 

San  Diego  Division 

271  Catalina  Boulevard 

San  Diego,  California  9^152 

Attn:  Code  d606 

Commander 

Naval  Undersea  Research  and  Development  Center 

San  Diego  Division 

271  Catalina  Boulevard 

San  Diego,  California  92152 

Attn : Library 

Commander 

Naval  Undersea  Research  and  Development  Center 
San  Diego  Division 
271  Catalina  Boulevard 
San  Diego,  California  92152 
Attn:  Bryson  Pennoyer 

Code  D554 

Commanding  Officer 

U.  S.  Fleet  ASW  School 

San  Diego,  California  921^7 

Commanding  Officer  and  Director 

U.  S.  Navy  Underwater  Sound  Laboratory 

Fort  Trumbull 

New  London,  Connecticut  O652O 
Attn:  Code  2242 

Director 

U.  S.  Naval  Research  Laboratory 
Department  of  the  Navy 
Anac ostia 

Washington,  D.  C.  20590 


i 


118 


June  1969 


DISTRIBUTION  LIST  FOR 
ARL-TM-69-11 

CONTRACT  NOOO24-69-C-II29  (CONT'D) 


Copy  No. 


Ik 

Commander 

U.  S.  Naval  Ordnance  Laboratory 
White  Oak 

Silver  Spring,  Maryland  2091O 

15 

Commanding  Officer 

U.  S.  Naval  Ship  Research  and  Development 
Laboratory'- 

Panama  City,  Florida  }2kOf. 

16 

Commanding  Officer 

U.  S.  Naval  Air  Development  Center 

Johns ville 

Warminster,  Pennsylvania  18974- 

17 

Commander 

Naval  Undersea  Research  and  Development  Center 
3202  East  Foothill  Boulevard 
Pasadena,  California  9HC7 

18 

Naval  Ship  Research  and  Development  Center 
Department  of  the  Na-vy 
Washington,  D.  C.  200C7 
Attn:  T,  Kooij 

19 

Office  of  the  Director  of  Defense 
Research  and  Engineering 
Room  3CI28,  The  Pentagon 
Washington,  D.  C,  20301 

20 

U.  S.  Na-val  Postgraduate  School 
Monterey,  California  93940 

21 

Commander 

Destroyer  Development  Group  TWO 
U.  S.  Atlantic  Fleet 
C/O  Fleet  Post  Office 
New  York,  New  York  09501 

119 


June  1969 


DISTRIBUTION  LIST  FOR 
ARL-TM-69-11 

CONTRACT  N00021)-  69-  C- 1129  ( CONT ' D ) 


25 

26 

27  - 51 

52 

53 

34  - 52 


Commander 

Submarine  Developnient  Group  TWO 
Box  70,  Naval  Submarine  Base 
New  London 

Groton,  Connecticut  065*4-0 

National  Academy  of  Sciences 
2101  Constitution  Avenue,  N.W. 
Washington,  D.  C.  20*4-18 
Attn:  Committee  on  Undersea  Warfare 

Office  of  Naval  Research 
Resident  Representative 
The  University  of  Texas  at  Austin 
2507  Main  Building 
Austin,  Texas  78712 

Signal  Physics  Division,  ARL/UT 

Computer  Science  Division,  ARL/UT 

J.  F,  Hoffman,  ARL/UT 

S.  P.  Pitt,  ARL/UT 

Library,  ARL/UT 

ARL  Reserve,  ARL/UT 


A 

I •' 


t - • 
?■ 


